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ABSTRACT 


No  matter  what  method  Is  employed  In  the  programming  of  a 
transducer  array,  in  practice,  the  performance  predicted  by  the  model 
will  not  be  achieved  due  to  error  in  construction  and  operation  of  the 
array.  Furthermore,  the  degree  to  which  a given  amount  of  error 
degrades  an  array  design  (i.e.,  design  sensitivity)  will  differ  greatly 
from  one  design  to  the  next.  A measure  of  sensitivity  is  proposed, 
based  on  the  maximum  impact  which  a given  amount  of  error  can  produce  on 
a particular  design.  The  adoption  of  this  measure  allows  the 
sensitivity  problem  to  be  solved  through  the  use  of  mathematical 
programming  techniques.  The  best  of  these  techniques  for  the  general 
array  design  case  is  shown  to  be  Non-linear  Goal  Programming.  The 
sensitivity  values  for  five  planar  symmetric  array  designs  are 
determined  and  compared. 

The  results  of  this  comparison  Indicate  the  potential  of  the 
proposed  measure  and  solution  technique  for  the  development  of  a fast 
and  Inexpensive  method  for  array  design  sensitivity  analysis. 
Particularly,  the  resulting  sensitivity  values  are  free  from  the  a 
priori  assumptions  and  high  costs  inherent  in  simulation.  Furthermore, 
the  method  is  based  on  a general  framework  which  is  applicable  to  the 
analysis  of  several  array  design  factors  including  size,  geometry, 
pattern,  and  component  reliability.  The  method  may  also  be  extended  to 
multiobjective  design  problems  outside  the  area  of  array  design. 
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CHAPTER  I 

( 


INTRODUCTION 


1.1  Purpose  of  the  Study 

Mathematical  programming  techniques  appearing  during  the  past  10 
years  have  provided  a means  for  designing  complex  array  patterns.  Prom 
1946  until  the  advent  of  these  techniques,  arrays  were  almost 

universally  programmed  using  the  Dolph  Tschebbyschef f method,  which 
Involved  setting  the  directional  response  function  equal  to  a 
Tschebbyschef f polynomial  of  like  degree. [1,9]  This  method  applied  only 
to  linear  arrays  In  the  exact  sense.  However,  by  the  use  of  the  product 

theorem,  patterns  produced  by  planar  arrays  could  be  approximated.  The 

Dolph-Tschebbyschef f method  allowed  main  beam  width  and  side  lobe  height 
to  be  specified.  In  contrast,  it  is  now  possible  to  achieve  main  beam 
shape  and  place  nulls,  as  well  as  achieve  other  objectives.  A 

concomitant  to  the  Increased  demands  on  array  performance  is  a decreased 
tolerance  of  the  array  pattern  to  small  deviations  in  the  values  of 
certain  key  parameters.  Thus,  the  potential  Impact  of  error  has  become 
of  increased  concern  to  the  designer.  This  concern  provides  the 
motivation  for  the  present  study,  the  purpose  of  which  is  to  find  a 
practical  means  of  assessing  array  pattern  sensitivity. 

This  process  entails  the  choice  of  some  measure  of  the  Impact  of 


2 


error  on  an  array  pattern  to  serve  as  an  indicator  of  sensitivity,  and 
concurrently,  the  development  of  a method  by  vhich  this  measure  can  be 
evaluated.  The  choice  reflects  a compromise  between  attaining  as  much 
information  regarding  sensitivity  as  possible,  and  keeping  the  time  and 
expense  Involved  in  its  acquisition  within  reasonable  limits. 

The  design  of  this  study  encompasses  five  main  objectives.  These 

are: 

1,  To  define  a consistent  measure  of  array  pattern 
sensitivity,  which  may  be  applied  to  all  classes  of  array,  and 
to  any  type  of  array  design, 

2,  To  develop  or,  more  correctly,  to  adapt  an  efficient 
method  to  determine  array  pattern  sensitivity,  based  on  the 
above  measure, 

3,  To  discuss  and  compare  alternate  means  for  measuring  and 
determining  sensitivity, 

4,  To  apply  the  chosen  method  to  a representative  problem. 
Specifically,  on  an  array  of  given  size  and  geometry,  the 
sensitivities  of  five  different  array  patterns,  displaying  a 
wide  range  of  possible  design  features,  will  be  determined  and 
compared, 

5,  To  discuss  the  extension  of  the  method  to  the  general 
array  design  sensitivity  analysis  problem. 

1,2  Organization  of  the  Paper 

The  paper  is  organized  as  follows.  Chapter  Two  includes  an 
introduction  to  transducer  arrays  and  array  patterns,  pertinent 
terminology,  coordinate  systems,  and  conventional  representations. 
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Methods  of  array  pattern  design  using  mathematical  programming 
techniques  are  summarized,  and  error  sources  are  discussed,  as  well  as 
their  resolution  into  variables  of  the  directivity  function* 

In  Chapter  Three,  alternate  methods  of  measuring  sensitivity  are 
discussed.  These  are  the  expected  value  and  the  maximum  value  of  change 
in  the  array  pattern.  Then,  methods  of  finding  these  values  are 
compared,  including  simulation  and  a set  of  mathematical  optimization 
techniques. 

MAXIM,  a program  which  measures  the  maximum  impact  of  error  on  an 
array  pattern  is  presented,  and  the  use  of  MAXIM  to  develop  an  index  of 
sensitivity  is  discussed. 

Chapter  Four  is  a description  of  an  experimental  application  of 
MAXIM.  Sensitivity  values  are  generated  for  five  directivity  patterns. 
The  results  are  presented  in  graphic  and  tabular  form. 

These  results  are  discussed  in  Chapter  Five,  and  extensions  and 
alternate  applications  of  the  method  are  suggested. 
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CHAPTER  II 


( 


ARRAY  DESIGN  AND  ERROR 


f 


2.1  Transducer  Arrays 

A transducer  is  any  device  which  converts  a signal  from  one 
transmission  medium  into  a signal  in  another.  When  operating  as  a 
reciever  the  transducer  converts  a signal  incident  to  its  surface  into 
an  electric  signal,  and  converts  an  electric  signal  into  a radiated 
signal  when  operating  as  a transmitter. 

A transducer  array  consists  of  a set  of  transducers,  termed 
elements,  which  are  placed  according  to  some  fixed  geometrical 
arrangement,  and  are  electrically  interconnected.  Arrays  fall  into 
three  categories , depending  on  the  number  of  physical  dimensions 
necessary  to  acomodate  the  element  positions.  These  are:  linear, 
planar,  and  conformal,  corresponding  to  one,  two,  and  three  dimensions, 
respectively. 

2.2  Array  Patterns 

One  of  the  major  objectives  Involved  In  acoustic  array  design  Is 
the  achievement  of  a desired  array  pattern,  or  equivalently,  directivity 
pattern.  If  the  operating  array  Is  viewed  from  an  aspect  angle, 
expressed  as  a solid  angle,  0,<),  where  the  0 component  Is  measured  from 
the  positive  Z axis,  and  the  <)  component  from  the  positive  X axis  (see 
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fig.2.1),  the  power  radiated  along  that  angle  (or  the  sensitivity  of  the 
array  to  a signal  Incident  to  It  from  that  angle)  Is  called  the 
directional  response  of  the  array  In  the  0,$  direction.  The  directional 
response  Is  measured  in  dB.  Letting  0 and  (}  vary  to  define  some  region 
in  space  (often  a hemisphere),  the  directional  response  at  all  angles 
contained  In  that  region  Is  called  an  array  pattern.  The  parameters 
which  determine  the  array  pattern  are: 

1.  The  number  and  placement  of  elements. 

2.  The  wavelength  at  which  the  system  operates. 

3.  The  relative  amplitude  and  phase  shading  coefficients 
of  the  excitation  currents  to  the  Individual  elements. 

Figure  2.2  Is  a contour  plot.  It  represents  a typical  array  pattern 

defined  on  0 < 0 < 180  and  0 < <1  < 180.  This  pattern  Is  obtained  from  a 

10  by  10  planar  symmetric  array  with  unshaded  element  amplitude  and 

phase  settings.  The  response  of  the  array  Is  given  on  a five-degree 

grid  In  negative  decibels.  Thus  -0  dB  represents  the  peak  power 

radiated  In  any  direction,  and  one  half  of  the  total  gain  (radiated 

power)  of  the  array  Is  in  the  region  -3  dB  or  less  off  peak.  In  figure 

2.2  the  center  of  this  region,  called  the  main  beam,  or  main  lobe.  Is  at 

0 - 90,  <)  • 90.  In  this  case,  the  axis  or  boreslght  of  the  array  falls 

along  the  positive  Y axis.  Any  other  lobe  In  the  array  pattern  with  a 

relatively  high  gain  , for  Instance  at  0 - 60,  <>  - 90,  is  considered  to 

be  noise  and  Is  called  a side  lobe.  Likewise,  any  region  where  the 

gain  falls  below  some  prespecified  level,  say  -60  dB , is  termed  a null, 

d 

e.g.,9  ■ 15,  0 - 40.  Figure  2«3  is  a simulated  3 plot  of  the  same 

pattern. 

An  array  pattern  may  be  defined  In  terms  of  three  categories  of 
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design  feature: 

1.  The  size  and  shape  of  the  main  beam. 

2.  The  "height"  of  the  side  lobes. 

3.  The  placement  of  nulls. 


2*3  Array  Design  Through  Optimization  Techniques 

During  the  past  ten  years,  problems  In  array  pattern  design  have 
been  successfully  solved  using  mathematical  optimization  techniques. 
Among  the  methods  which  have  been  employed  are  Lagranglan  multipliers, 
linear  programming,  method  of  steepest  ascent,  SUMT,  and  noA-linear  goal 
programming.  Although  these  methods  differ  widely  with  regard  to 
solution  technique  (and  relative  success),  they  usually  share  an 
underlying  view  of  the  array  design  problem  as  an  optimization  model  In 
m objectives,  and  2n  variables.  Thus  an  array  of  n independent 
elements,  or  sets  of  elements,  is  required  to  meet  a set  of  specified 
directional  responses,  at  m aspect  angles  (0,  0 points),  thereby 
defining  an  array  pattern.  In  general.  It  Is  desired  to  maximize  an 
objective  function 

Z(fl(al a2n) ,fm(al a2n}’  f'l“**’f'm) 

where  f^fa^ ,. . . ,a2n)  18  thc  achieved  directional  response  at  0^,0^, 
a^,...a2n  are  the  amplitude  and  phase  shading  coefficients  of  the 
elements  of  the  array,  and  f'^  Is  the  specified  response  at  9^,0^  The 
function  Z will  take  on  Increasingly  greater  values  as 


fjUj.***  i*2n)  ->  f#t»  1-1,... ,« 
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Due  to  a somewhat  unfortunate  overlap  in  terminology,  the  set  of  values 
which  Z takes  on  as  ap...,a2n  vary  continuously,  Is  called  a response 
surface.  The  term  "response  surface"  has  to  do  with  the  response  of 
the  objective  function  and  does  not  refer  directly  to  the  directional 
response  of  the  array  pattern.  A response  surface  In  2 rather  than  2n 
dimensions  is  shown  in  figure  2. A. The  point  (a'^.a'^)  Is  called  a global 
optimum.  That  Is,  there  is  no  point  (a^a^)  in  the  domain  of  Z such 
that  Z(a1,a2)  > Z(a'1,a'2>.  The  point  (a"  ,a"2)  is  a local  optimum, 
since,  although  the  conditions  for  global  optimality  do  not  hold,  there 
exists  a region  given  by: 

I < di»  d±  > 0 

such  that  Z(a^,a2>  _<  Z(a"^,a"2)  holds  for  all  a^,a2  within  the  region. 

Of  the  optimization  methods  listed  above,  Lagrangian  multipliers 
and  linear  programming  are  called  indirect  methods,  since  their  solution 
techniques  involve  the  simultaneous  solution  of  a set  of  equations.  The 
method  of  steepest  ascent,  the  sequential  unconstrained  minimization 
technique,  and  non-linear  goal  programming,  are  direct  search  techniques 
relying  on  evaluation  of  response  surface  elevation,  and  some  search 
movement  logic  to  find  an  optimum. 

Recently  [10,15],  Non-linear  Goal  Programming  (NLGP ) , which  is 
based  on  the  modified  Hooke  and  Jeeves  pattern  search  discussed  in 
appendix  B,  has  been  shown  to  provide  good,  robust  solutions  to  the 
array  design  problem.  The  success  of  NLGP  is  due  in  part  to  the 
simplicity  and  directness  of  its  search  technique.  All  necessary 
Information  for  search  movement  decisions  is  available  directly  from  the 
objective  function.  Another  factor  in  the  success  of  NLGP  derives  from 
the  goal  programming  problem  formulation,  which  accommodates  multiple 
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possibly  conflicting,  or  even  lncommeasurable  objectives.  This 
particular  feature  of  NLGP  is  especialy  applicable  to  the  array  design 
problem,  since  the  achievement  of  different  array  pattern  features 
usually  conflict.  For  example,  it  is  difficult  to  achieve  a narrow  main 
beam,  and  at  the  same  time  maintain  low  side  lobe  levels. 

2.4  Array  Sensitivity  and  Error 

Sensitivity  analysis  is  the  measurement  of  system  response  to 
deviations  in  the  values  of  decision  variables  from  those  assumed  in  the 
optimal  solution.  The  technique  by  which  sensitivity  is  found  is  not 
dependent  on  the  method  originally  used  to  obtain  the  solution. 

If,  as  is  often  the  case,  the  size  , shape,  and  operating 
wavelength  of  the  array  have  been  specified  in  advance,  the  array  design 
solution  consists  of  a vector  of  element  amplitude  and  phase  shading 
coefficients.  It  happens  that  the  values  of  these  variables,  specified 
in  the  solution  are  not  precisely  attained  in  practice.  The  reason  for 
this  may  be  traced  to  several  sources  of  error,  inherent  in  the 
construction  and  operation  of  a transducer  array.  The  following  list  of 
error  sources  is  not  intended  to  be  comprehensive,  but  to  provide  a 
representative  sample. 

Error  sources  are  usually  divided  into  two  classes , depending  on 
whether  their  effects  on  the  array  pattern  are  random  or  systemmatic. 

The  most  important  error  source  in  the  random  category  is  control 
error  in  the  excitation  signals  recleved  by  the  individual  elements. 
This  error  source  is  analogous  to  instrument  drift,  and  will  vary  with 
time.  Random  construction  errors  in  the  manufacture  and  placement  of 
elements,  in  the  manufacture  of  phase  shifters,  or  in  the  length  and 
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fitting  of  feed  lines  will  have  a fixed  but  unpredictable  effect  on  a 
given  array. 

Chief  among  the  systemmatlc  errors  is  a phenomenon  called  mutual 
coupling.  The  signal  emitted  from  any  one  element  will  induce  a 
sympathetic  excitation  in  every  other  element  in  the  array.  The 
amplitude  and  phase  of  this  component  of  a given  element's  radiation 
depend  on  its  position  in  the  array,  and  the  amplitude  and  phase 
setting  of  each  other  element.  Systematic  phase  error  may  also  be 
introduced  by  sensitivity  of  phase  shifters,  transmission  lines,  and 
other  array  components  to  frequency  variations.  Another  source  of 
phase  error  is  the  approximation  Inherent  in  discrete  or  digital  phase 
shifters. 

The  effects  of  these  and  any  other  error  source  may  be 
represented  as  a change  in  element  amplitude  and  phase  - shading 
coefficients.  However,  the  error  Introduced  by  the  mathematical  model 
of  directional  response,  which  assumes  elements  with  point  apertures, 
and  isotropic  radiation  patterns  is  not  treated  here.  This  error  source 
may  be  resolved  through  refinements  in  the  model  if  necessary. 
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CHAPTER  III 


THE  MEASUREMENT  OF  SENSITIVITY 


3.1  Sensitivity  by  Feature 

In  the  array  design  problem  the  analyst  is  forced  to  consider  the 
entire  pattern  at  once,  which  results  in  a multiobjective  problem 
formulation.  This  is  because,  as  has  been  previously  noted,  the 
achievement  of  different  array  pattern  features  often  conflict.  While  a 
pattern  satisfying  all  design  objectives  would  nominally  be  optimal,  it 
probably  won't  be  feasible.  Optimality  is  then  assigned  to  the  "best" 
trade  off  which  is  both  acceptable  and  feasible.  The  relative 
importance  of  different  features  may  be  reflected  either  in  a set  of 
weighting  factors  or  in  a preemptive  priority  structure.  It  is 
precisely  this  type  of  multiobjective  problem  which  goal  programming  was 
designed  to  solve. 

In  the  case  where  array  pattern  sensitivity  is  to  be  determined, 
however,  it  is  more  profitable  to  examine  individual  features  of  an 
array  pattern  separately.  This  is  because  the  degradation  of  the 
different  features  will  conflict  in  exactly  the  same  manner  as  their 
achievement.  By  separating  a pattern  with,  say,  four  important  features 
into  four  different  sensitivity  problems,  the  question  of  deciding  the 
relative  "worst"  between  alternative  examples  of  impact  where  different 
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features  are  affected  is  circumvented.  With  this  separation,  four  new 
response  surfaces  may  be  created,  each  representing  the  achievement  of 
one  feature.  Although  the  solution  vector  lies  at  the  optimal  point  of 
the  array  pattern  response  surface,  it  does  not  follow  that  it  lies  at 
the  optimal  point  of  any  of  the  individual  feature  response  surfaces. 


3.2  The  Error  Region  and  Error  Surface 

Consider  the  solution  vector  of  an  array  design  problem.  Each 
component  of  this  vector  may  be  said  to  have  an  interval  of  error  about 
it,  representing  the  potential  cumulative  effects  of  all  error  sources 
which  can  be  resolved  into  error  in  amplitude  and  phase  - shading 
coefficients.  The  true  value  of  the  elements  will  lie  somewhere  along 
these  Intervals,  which  are  expressed  as  plus  or  minus  some  percentage  of 
the  range  over  which  the  amplitude  or  phase  can  vary.  Amplitudes  are 
conventionally  normalized  to  dimensionless  values  between  0 and  1,  phase 
to  values  between  0 and  360  degrees.  For  an  n-element  array,  these 
error  Intervals  will  form  a hypercube  in  2n  space,  with  the  solution 
vector  at  its  center.  This  hypercube  will  be  called  the  error  region. 
A 2-dimensional  error  region  is  shown  in  figure  3«1. 

If  the  error  region  is  superimposed  on  a feature  response  surface, 
it  will  contain  all  values  which  that  objective  function  which  indicates 
the  achievement  of  the  desired  feature  will  take  on  due  to  error.  This 
part  of  the  response  surface  will  be  called  the  error  surface. 

Perfect  information  about  the  sensitivity  of  a feature  would 
consist  of  two  types  of  data.  The  first  of  these  would  be  a complete 
map  of  the  error  surface.  This  is  equivalent  to  complete  enumeration  of 
the  response  function  value  over  2n  continuous  Intervals.  It  is  easily 
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seen  that  even  a low-resolution  approximation  would  require  that  an 
astronomical  number  of  points  be  evaluated  as  the  array  becomes  large. 
For  example,  the  array  pattern  of  a 100  -element  array  depends  on  200 
decision  variables.  If  the  error  interval  of  each  variable  is 
approximated  by  50  equally  spaced  points,  forming  a 200  -dimensional 
lattice,  total  enumeration  of  the  error  surface  response  at  the  lattice 
points  requires  the  array  pattern  to  be  evaluated  10,000  times.  It 
should  be  noted  as  well,  that  as  the  array  size  increases,  so  does  the 
cost  of  each  evaluation.  The  second  set  of  data  would  be  composed  of 
detailed  knowlege  of  the  effects  of  error  sources.  This  would  Include 
such  information  as  analysis  of  component  interaction  in  various  array 
constructions,  probability  distributions,  the  different  random  error 
types,  and  so  forth.  This  kind  of  data  is  not  available  in  practice. 

Even  though  perfect  sensitivity  information  is  impractical  or 
impossible  to  obtain,  it  is  still  desirable  to  find  some  characteristic 
of  the  error  surface  which  will  allow  an  assessment  of  sensitivity  to  be 
made.  Candidates  include  the  average  elevation  and  the  range  (maximum 
elevation  minus  minimum  elevation)  of  the  error  surface.  The 
corresponding  measures  would  be,  respectively,  the  expected  value  of  the 
feature's  performance,  and  the  worst  value  of  the  feature's  performance 
in  the  error  region. 

3.3  Simulation 

An  estimate  of  the  mean  and  standard  deviation  of  the  error  surface 
elevation  can  be  obtained  by  sampling  points  within  the  error  region, 
and  evaluating  the  objective  function  at  these  points.  This  amounts  to 
a simulation  of  error  impact.  The  sample  will  consist  of  a set  of 
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points,  constrained  to  lie  within  the  error  region,  and  randomly 
generated  according  to  some  distribution.  For  this  purpose  the 
distribution  will  usually  be  either  normal  about  the  solution  vector, 
with  a specified  standard  deviation,  or  uniform  across  the  error  region. 
Since  the  error  sources  discussed  above  are.  In  fact,  partially 
systematic,  and  their  random  components  partially  fixed  with  regard  to 
effect  on  a particular  array,  the  assumptions  Inherent  In  these  types  of 
Monte  Carlo  simulations  may  seriously  bias  the  analysis. 

A second  problem  with  simulation  is  the  required  sample  size,  and 
consequent  expense.  The  number  of  degrees  of  freedom  In  the  directivity 
function  is  twice  the  number  of  Independent  elements  in  the  array. 
Reliable  sample  sizes  tend  to  become  quite  large.  It  Is  noteworthy  that 
in  the  100-element  case  mentioned  earlier  (not  a particularly  large 
array),  even  a ten  per  cent  sample  amounts  to  1000  points. 

3*4  Optimization  Models 

Since  simulation  and  enumeration  are  the  only  means  of  arriving  at 
the  expected  value  of  array  pattern  feature  sensitivity,  and  both  prove 
to  be  impractical  in  light  of  present  computer  speeds  and  capacities,  it 
remains  to  examine  the  error  surface  range.  The  maximum  elevation  on 
the  error  surface  is  assumed  to  be  the  solution  vector  for  the  array 
pattern.  Quite  a few  optimization  methods  may  be  used  to  find  the 
minimum  elevation.  These  turn  out  to  be  the  same  techniques  which  have 
previously  been  brought  to  bear  on  the  array  design  problem.  It  should 
not  be  surprising,  considering  the  relationship  between  the  two 
problems,  that  the  relative  success  of  a technique  In  solving  one  is  a 
good  predictor  of  performance  on  the  other. 
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Representatives  of  three  classes  of  optimization  models  will  be 
discussed.  These  are: 

1.  Calculus-based  methods. 

2.  Linear  Programming. 

3*  Direct  search  techniques. 

In  order  to  restrict  an  Indirect  method,  l.e.  one  based  on  the 
solution  of  a set  of  equations,  to  the  error  region,  two  Inequality 
constraints  must  be  added  to  the  problem  for  each  decision  variable.  A 
typical  calculus-based  technique  for  minimization  subject  to  constraints 
is  the  use  of  Kuhn  Tucker  conditions.  This  method  uses  Lagrangian 
multipliers  to  create  a set  of  equations  which  are  solved 
simultaneously. 

In  the  case  of  the  general  (non-symmetrlcal)  n-element  array  these 
consist  of: 

1.  2n  Lagranglans  of  the  directivity  function  set  equal  to  0. 

2.  4n  Inequality  constraints. 

3.  4n  slack  constraints,  bounding  the  values  which  the 
multipliers  may  take  on. 

To  find  the  mlnumum  on  a region,  different  combinations  of  the 

multipliers  are  set  to  0,  forcing  the  method  to  consider  stationary 

points  only  along  one  or  more  boundaries  of  the  error  region.  No 

multipliers  set  to  0 corresponds  to  the  Interior  of  the  error  region. 

This  process  Is  comblnatorlcally  explosive.  At  the  worst  case,  up  to 
An 

2 separate  problems,  each  non-linear  In  2n  variables  and  lOn 
constraints,  may  have  to  be  solved.  Solving  even  one  of  these  problems 
however,  can  be  a somewhat  monumental  undertaking. 

Two  formulations  which  linearize  the  directivity  function  have  been 
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developed.  McMahon  et  al. [18]  were  the  first  to  translate  the  array 
pattern  design  problem  Into  a linear  programming  problem.  This  approach 
will  solve  for  both  amplitude  and  phase  coefficients,  but  in  order  to  do 
so,  is  forced  to  separate  the  real  and  imaginary  components  of  the 
pattern.  Because  of  this,  it  Is  impossible  to  reduce  the  problem  size 
by  taking  symmetry  Into  account. 

Unlike  calculus-  based  techniques,  which  are  applied  to  the  error 
surface,  linear  programming  evaluates  the  directivity  pattern.  Thus  a 
set  of  9,  <)  points  must  be  specified  to  define  a feature  which 
corresponds  to  some  worst  case.  At  each  9,  0 point,  four  Inequalities 
are  generated,  two  each  for  the  real  and  Imaginary  parts  of  the  pattern. 
Two  synthetic  decision  variables  are  developed  for  each  element  In  the 
array,  but  as  these  are  unrestricted,  their  number  Is  functionally 
doubled.  Also,  for  each  element,  four  Inequality  constraints  are  added 
In  order  to  define  the  error  region. 

Even  before  considering  the  number  of  9,  0 points  required  to 
define  the  worst  case  for  a feature,  there  will  be  a problem  In  4n 
variables,  and  4n  constraints.  While  this  problem  Is  orders  of 
magnitude  smaller  than  that  generated  by  Kuhn  Tucker  conditions.  Its 
solution  Is  still  an  expensive  project.  There  will  also  be  a certain 
amount  of  unavoidable  round  off  error  due  to  the  subtraction  and 
division  required  at  each  pivot. 

In  a later  linear  formulation  , Wilson  [20]  created  a linear 
programming  problem  which  deals  only  with  the  real  component  of  the 
pattern.  Symmetry  conditions  may  be  Incorporated  Into  the  directivity 
function  In  this  model,  and  the  decision  variables  are  restricted  to 
non-negative  values.  However,  this  drastic  reduction  In  problem  size  is 
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accomplished  at  the  expease  of  the  abilty  to  solve  for  phase  settings.  i 

« 

4 

Wilson's  model  assumes  that  all  elements  are  excited  in  phase,  and 

I 

solves  only  for  the  amplitude  coefficients.  Although  this  technique  j 

could  be  easily  adapted  to  a sensitivity  analysis  of  the  restricted 
(planar  symmetric  and  unphased  array)  class  of  problem  which  it  was 
intended  to  solve,  it  is  not  applicable  to  the  general  case. 

Direct  search  techniques  have  been  the  most  effective  means  of  < 

solving  the  general  array  design  problem.  The  sequential  unconstrained 
minimization  technique,  (SUM!)  [11,  12],  gradient  search  (method  of 

steepest  ascent)  , and  Non-linear  Goal  Programming  [4]  have  all  been 

utilized  to  search  for  an  optimum  directivity  pattern.  Summaries  and  ! 

comparisons  of  these  techniques  have  been  made  by  Ignlzlo  [15]  and  by 

| 

Draus  [10].  When  applied  to  Identical  pattern  synthesis  problems,  NLGP  \ 

was  found  to  equal  or  exceed  (sometimes  markedly)  the  performance  of  i 

the  other  two  methods,  while  at  the  same  time  displaying  a clear 

. 

superiority  in  computational  efficiency.  These  results,  plus  the 

multiobjective  capacity  of  the  goal  programming  formulation  make  NLGP  < 

the  most  attractive  direct  search  technique,  and  by  extension,  the  best 
optimization  model  to  adapt  to  the  general  directivity  pattern 

' 

sensitivity  problem. 


3.5  MAXIM/PS 

NLGP /PS  is  the  goal  programming  code  specifically  designed  for 
planar  symmetric  array  pattern  design.  This  code  was  written  by  James 
P.  Ignlzlo  and  Susan  M.  Draus  at  the  Industrial  Engineering 
Department,  and  Applied  Research  Laboratory,  The  Pennsylvania  State 
University.  Its  accuracy  .ind  robustness  have  been  extensively  tested. 
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Only  minor  revisions  were  necessary  in  order  to  adapt  NLGP  to  search  for 
the  minimum  elevation  of  the  error  surface.  These  changes  involve 
parameters  of  the  pattern  search,  which  is  described  in  detail  in 
appendix  B. 

1.  The  initial  step  sizes  were  reduced  to  compensate  for  the 
constrained  region  in  which  the  search  takes  place. 

2.  Logic  was  added  to  the  evaluation  of  provisional  base 
points  such  that  any  point  outside  the  error  region  registers 
no  improvement.  This  effectively  constrains  the  search  to  the 
error  region. 

3.  Changes  were  made  in  the  functional  relationships  which 
define  the  goals,  that  is,  the  response  desired  at  given  0,  <> 
points.  In  NLGP  the  achieved  levels  of  response  were  set 
equal  to  the  desired  levels  in  the  main  region,  and  less  than 
or  equal  in  the  side  region.  In  the  revised  version,  response 
in  the  main  region  is  set  less  than  or  equal  to  the  desired 
level,  and  in  the  side  region,  greater  than  or  equal  to  the 
desired  level. 

4.  The  solution  vector  of  the  design  problem  is  used  as  the 
initial  base  point  for  the  pattern  search. 

The  goal  programming  formulation,  the  directional  response 
formulation,  and  the  pattern  search  algorithm  are  identical  in  both 
codes.  The  revised  version  is  called  MAXIM/PS,  as  it  measures  the 
MAXlmum  IMpact  of  an  error  range  on  a feature. 

Direct  search  techniques  do  not  guarantee  a global  optimal  solution 
on  a multimodal  response  surface,  and  the  Hooke  and  Jeeves  pattern 
search  is  no  exception.  This  is  less  of  a problem  in  the  array  design 


case,  because  only  a satisfactory  solution  Is  required.  The  question  of 
whether  or  not  the  global  minimum  of  an  error  surface  has  been  found  Is 
of  greater  concern  , since  a difference  of  several  dB  of  sensitivity  may 
be  involved.  In  practice,  however,  error  Intervals  are  very  small  (on 
the  order  of  + 2 or  3 percent  of  the  decision  variable  range),  so  that 
the  chance  of  multimodality  on  any  error  interval,  and  hence  on  the 
error  surface.  Is  significantly  reduced.  If  an  extra  measure  of 
confidence  is  required,  it  is  possible  to  make  multiple  search  runs, 
commencing  from  different  starting  points. 

3.6  The  Index  of  Sensitivity 

As  noted  previously,  in  the  determination  of  pattern  sensitivity 
Important  features  should  be  considered  separately.  Thus,  the  9,  0 
points  corresponding  to  a high  side  lobe  might  be  chosen  and  assigned 
an  aspired  value  of  -1  dB.  MAXIM  would  then  search  within  the  error 
region  for  the  point  at  which  this  side  lobe  was  driven  as  high  as 
possible.  A measure  on  the  9,  4 points  used  to  define  the  side  lobe, 

V 

say  the  highest  point  or  the  mean  height,  could  then  be  compared 
between  the  solution  to  the  design  problem  and  the  maximum  error  Impact 
solution.  If  the  directional  response  at  these  Q , 0 points  had  risen  by 
an  average  of  15  dB,  or  if  the  maximum  height  of  the  feature  had  risen 
by  15  dB,  the  feature  would  be  given  a sensitivity  value  of  15  dB.  The 
same  procedure  would  then  be  repeated  for  all  relevant  features  of  the 
pattern.  These  would  Include  other  side  lobes,  nulls,  and  the  main 
beam.  The  main  beam  9,  0 points  would  be  given  an  aspired  value  of  -60 
dB  to  bring  the  response  down  as  far  as  possible. 

For  a given  error  range,  the  result  of  this  approach  Is  a set  of 
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sensitivity  values  associated  with  the  Important  features  of  the  array 
pattern.  These  are  to  be  compared  with  minimum  acceptable  sensitivity 
values  for  these  features.  The  maximum  amount,  in  dB , by  which  the 
acceptable  sensitivity  value  for  any  feature  of  the  pattern  is  violated 
is  the  I OS , or  index  of  sensitivity  assigned  t»  the  pattern.  Note  that 
the  IOS  is  specific  to  a given  error  range. 


CHAPTER  IV 


( 


APPLICATION:  THE  USE  OF  MAXIM  TO  GENERATE  AN  IOS 


4* 1 The  Experiment 

In  order  to  evaluate  the  utility  of  the  IOS  measure  of  sensitivity, 
it  was  decided  to  test  and  compare  a series  of  directivity  patterns  on  a 
single  array.  The  directivity  patterns  were  chosen  so  as  to  reflect  a 
broad  sampling  of  objectives.  Synthesis  took  place  on  a 10  x 10  planar 
symmetric  array.  The  choice  of  array  size  and  geometry  was  based  on  two 
considerations.  The  number  of  elements  was  to  be  large  enough  to 
produce  patterns  which  were  somewhat  difficult,  and  hence  interesting 
from  the  standpoint  of  sensitivity.  At  the  same  time,  the  number  of 
decision  variables  involved  was  to  be  kept  down  to  a reasonable  level, 
allowing  the  analysis  to  be  accomplished  at  a small  to  moderate  cost  in 
computation  time. 

4.2  Array  Description 

The  10  x 10  planar  symmetric  array  is  square  and  consists  of  100 
f equally  spaced  elements.  The  elements  are  considered  to  lie  in  the  X-Z 

plane,  25  elements  to  a quadrant,  with  symmetry  about  both  the  X and  Z 
axes.  The  elements  are  placed  and  numbered  as  shown  in  figure  4.1.  In 
the  planar  symmetric  type  of  array,  sets  of  four  symmetrically  placed 

r 

| 

T-:r  j ■m*-  ■ 

, . - .*$«*(-* 


27 


elements  recleve  identical  excitation  currents.  This  fact  decreases 
both  the  amount  of  hardware  necessary  for  array  construction,  and  the 
difficulty  in  programming  the  array.  In  the  10  x 10  case,  it  is 
necessary  to  consider  only  25  sets  of  elements,  a reduction  from  200  to 
50  decision  variables. 

The  directivity  pattern  is  measured  over  a hemisphere  defined  by: 

0 < 9 < 180  and  0<  $ < 180 


See  figure  2*1.  Thus  boresight  falls  along  the  positive  Y axis. 
The  directivity  function  for  this  array  is  given  by: 


25  2 25 

10  Log  [16  [ (E  A.  cos  a cos  b cos  c ) + (E  A sin  a cos  b. 

i*l  1 1 1 1 i-1  1 


cos  c^  ]] 


where: 

a 4 ■ the  phase  shift  associated  with  element  1. 
b^^  - 2'n/XXjSin9cos0 

c^  ■ 2ir/Xzicose 

A^  - the  excitation  amplitude  of  element  1. 

X - the  operating  wavelength  of  the  array. 
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See  appendix  A for  the  derivation  of  this  result  from  the  general 
directivity  function. 

4.3  Experimental  Patterns 

Five  patterns  were  synthesized  using  NLGP.  This  process  is 
documented  in  references  10  and  15.  These  are: 

1.  20  by  30  degree  elliptical  main  beam  with  uniform 
response.  Side  lobes  _<  -27  dB  off  peak.  -61  dB  nulls  at 
(20,20),  (160,20),  (20,160),  (160,160). 

2.  30  degree  circular  main  beam  with  uniform  response.  Side 
lobes  £ -26  dB  off  peak. 

3.  30  degree  circular  main  beam  with  uniform  rasponaa.  Side 

lobes  £ -26  dB  off  peak.  Broad  nulls  _<  -47  dB  off  peak, 
centered  about  (40,35),  (40,165),  (140,165),  (140,35), 

(140,165). 

4.  30  degree  square  main  beam  with  uniform  rasponaa.  Slda 
lobes  -28  dB  off  peak. 

5.  30  by  70  degree  rectangular  main  beam  with  uniform 
response.  Side  lobes  < -27  dB  off  peak. 

In  the  discussion  to  follow,  the  experimental  patterns  will  be  rafarrad 
to  as  patterns  1 through  5.  These  five  patterns  are  illustrated  by  5 
degree  resolution  contour  plots,  and  2 degree  resolution  simulated  3d 
plots  in  figures  4.2  to  4.11.  Tables  4.12  to  4.15  show  the  element  set 
amplitude  and  phase-shading  coefficients  for  the  five  design  solutions. 
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Figure  4.2 

Pattern  1.  Design  Solution 
Simulated  3**  Plot 
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Figure  4.11 

Pattern  5.  Design  Solution 
Contour  Plot 
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TABLE  4.1 


Pattern  1.  Amplitude  and  Phase  settings. 


Design  Problem  Solution  Vector. 


Element  Number 


1 

1.00000 

2 

0.83026 

3 

0.55588 

4 

0-28236 

5 

0-14425 

6 

0.73256 

7 

0.55325 

8 

0.36320 

9 

0.17836 

10 

0-00016 

11 

0.47729 

12 

0.31261 

13 

0.16010 

14 

0.06544 

15 

0.03629 

16 

0.35335 

17 

0.24252 

18 

0.18708 

19 

0.00300 

20 

0.05679 

21 

0.19194 

22 

0.16094 

23 

0.09302 

24 

0.02788 

25 

0.01390 

45.08598 

50.47632 

61.58244 

81.43579 

125.26320 

29.22868 

40.92334 

46.62827 

86.20399 

-23.47665 

5.06743 

4.35339 

7.01000 

-72.95085 

-110.93900 

-27.58368 

-58.43025 

-57.70381 

-96.97165 

-95.49803 

-79.03387 

-56.59669 

-116.47880 

-106.64440 

-153.24800 
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TABLE  4.2 


Pattern  2.  Amplitude  and  Phase  settings. 
Design  Problem  Solution  Vector. 


Element  Number 


1.00000 

0.72492 

0.40101 

0.21527 

0.16095 

0.67942 

0.49048 

0.23734 

0.24548 

0.13429 

0.24719 

0-21164 

0.27099 

0-21564 

0.13460 

0.14507 

0.14949 

0.21673 

0.14453 

0.02817 

0.13854 

0.11909 

0.06781 

0.07849 

0.00075 


86.59134 

76.29446 

46.41832 

-16.98532 

-59.23582 

81.91975 

63.05440 

18.64299 

-33.18059 

-52.38556 

58.71460 

19.38586 

-8.95276 

-35.90080 

-62.06729 

-31.14095 

-28.63022 

-37.39592 

-54.42915 

-68.23637 

-65.68265 

-53.94441 

-47.95035 

-64.44748 

-78.08585 
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TABLE  4.3 

Pattern  3.  Amplitude  and  Phase  settings. 
Design  Problem  Solution  Vector. 


Element  Number 

Ai 

ai 

1 

1.00000 

86.58992 

2 

0.72632 

76.29305 

3 

0.34581 

51.41776 

4 

0.23573 

-21.98474 

5 

0.14923 

-43.99036 

6 

0.69419 

81.91833 

7 

0.47213 

63.05298 

8 

0.24763 

18.64278 

9 

0.23630 

-30.62059 

10 

0.15778 

-61.38362 

11 

0.24962 

58.71346 

12 

0.19205 

21.43350 

13 

0.26085 

-8.95269 

14 

0.20757 

-40.90024 

15 

0.10105 

-52.45972 

16 

0.14143 

-34.98846 

17 

0.14389 

-26.19022 

18 

0.20862 

-37.39572 

19 

0.13913 

-49.42932 

20 

0.05563 

-57.98967 

21 

0.11944 

-68.87962 

22 

0.11464 

-57.30469 

23 

0.06528 

-43.59050 

24 

0.07555 

-50.84010 

25 

0.00073 

-73.46278 
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TABLE  4.4 

Pattern  4.  Amplitude  and  Phase 

settings. 

Element  Number 

Design  Problem  Solution  Vector. 

Ai  “i 

1 

1.00000 

83.71129 

2 

0.58462 

74.89516 

3 

0.14241 

34.26509 

4 

0.19317 

-61.69458 

5 

0.10087 

-100.38940 

6 

0.55779 

84.17863 

7 

0.34828 

59.88213 

8 

0.23030 

6.18118 

9 

0.18049 

-23-08340 

10 

0.10228 

-77.36714 

11 

0.00410 

-174.76480 

12 

0.14119 

-0.87991 

13 

0.26177 

5.56305 

\ 

I 

14 

0.19908 

-13.63951 

15 

0.08626 

-12.35343 

• 

16 

0.15645 

-99.23828 

17 

0.10658 

-45.02731 

18 

0.16667 

-6.84196 

19 

0.14051 

1.34623 

1* 

20 

0.01905 

-10.38139 

' 

21 

0.08704 

-100.32360 

- j 

22 

0.03587 

-78.75728 

23 

0.03556 

-40.04144 

4 

24 

0.02226 

115.26540 

p 

25 

0.01850 

89.00087 
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TABLE  4.5 

Pattern  5.  Amplitude  and  Phase  settings. 
Design  Problem  Solution  Vector. 


Element  Number 

- Ai 

1 

0.66059 

102.85650 

2 

0.38303 

110.02400 

3 

0.01744 

105.41250 

4 

0.09856 

-85.83286 

5 

0.07641- 

-48.83804 

6 

1.00000 

41.89236 

7 

0.65056 

37.02228 

8 

0.17021 

38.25259 

9 

0.08036 

-93.32645 

10 

0.14957 

-147.13670 

11 

0.81610 

2.49893 

12 

0.47947 

1.34707 

13 

0.17213 

-20.41713 

14 

0.07275 

-167.79430 

15 

0.06579 

-165.56710 

16 

0.53865 

-75.06853 

17 

0.36734 

-77.70242 

18 

0.12291 

-93.52895 

19 

0.00422 

63.56360 

20 

0.08095 

68.92059 

21 

0.27639 

-101.25180 

22 

0.17685 

-107.32960 

23 

0.05861 

-94.70036 

24 

0.05308 

43.36661 

25 

0.03110 

53.67607 
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f For  each  pattern,  sensitivity  values  were  determined  for: 

1.  The  main  beam, 

2.  the  two  highest  side  lobes, 

' 3.  nulls  if  present, 

at  + 1 percent  and  + 2 percent  error  in  A^  and  & ^ values , and  at  -2 

percent , -1  percent , 1 percent  and  2 percent  error  in  1 . Because  of 

the  symmetry  inherent  in  the  pattern,  a side  lobe  or  null  in  one 
orthont  will  have  an  identical  counterpart  in  the  other  three.  For 
convenience,  a set  of  four  such  symmetrical  features  will  be  referenced 
only  in  the  orthont  defined  by: 

0 £ © £ 90  and  0 £ ()  £ 90 

The  process  by  which  the  sensitivity  values  were  obtained  is 
described  in  detail  for  pattern  1.  It  holds  for  the  other  four  patterns 
as  well. 
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4.4  The  MAXIM  Runs 

Four  major  features  were  identified.  These  are  the  main  beam,  the 
two  highest  side  lobes  , denoted  sl(l)  and  sl(2),  and  the  nulls.  These 
features  are  outlined  in  figure  4.2. 

The  sensitivity  of  each  of  these  features  to  error  in  A^  and  o was 
analyzed  separately,  using  MAXIM.  Thus  two  MAXIM  runs  were  made  for 
each  feature,  with  error  ranges  for  A^  and  set  first  at  1 percent, 
and  then  at  2 percent.  A^  and  vary  Independently,  and  the  capability 
to  model  this  is  built  into  MAXIM.  Setting  these  error  ranges  equal 
reflects  a desire  to  keep  the  analysis  as  simple  as  possible. 

Since  MAXIM  evaluates  the  directivity  pattern,  rather  than  the 
response  surface,  it  is  necessary  to  develop  a set  of  goals  which  will 


•W* 


drive  the  solution  to  the  minimum  elevation  point  on  the  error  surface. 
Directional  response  values  which  correspond  to  extreme  degradation  are 
therefore  specified  at  a number  of  9,  ( points  in  the  feature.  It  was 
felt  that  a 5-degree  resolution  was  sufficient  to  define  the  features 
examined  here.  Smaller,  irregularly  shaped  features  may  require  a finer 
resolution.  The  goals  are  of  the  form: 

Ge<f>  + ne<t>  " pe<j>  * Be<t> 


where: 


for  9,0  in  the  main  beam 


be 
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for  9,0  in  a side  lobe  or  null 

where  the  weighting  factors,  Wg^,  are  set  to  unity.  Indicating  no 
difference  In  priority  between  0,  0 points  In  the  same  feature.  This  is 
the  basic  goal  programming  formulation  in  one  priority  level.  There 
are,  in  fact,  three  priority  levels  in  the  problem,  but  two  of  these, 
excitation  symmetry  and  error  range,  are  Implicit  In  the  directivity 
function  and  pattern  search  logic  respectively. 

The  solution  vector  of  the  design  problem  was  used  as  the  initial 
base  point  for  the  pattern  search. 


4.5  Sensitivity  Values 

MAXIM  generates  a contour  plot,  displaying  the  worst  possible 
response  of  the  given  feature  on  a specified  error  region.  This  Is  then 
compared  with  the  contour  plot  of  the  design  solution  directivity 
pattern,  previously  generated  by  NLGP.  Two  sensitivity  values  are  taken 
for  each  MAXIM  contour  plot,  l.e.,  for  each  feature  at  each  error  range. 
These  are  the  average  difference  (AG),  and  the  difference  In  the  worst 
directional  response  values  for  any  9,  0 point  in  the  feature,  (AGmax). 

AG  is  given  by: 


1 (Ga*  • Ga*> 

e<p  -M si 


for  9,0  points  in  a side  lobe  or  null.  Or: 


Z (G 
60  - 


90 


I 


for  9,0  points  in  the  main  beam 


( 


Where 

G » Che  response  at  9,(1)  for  the  design  solution* 

0<j> 

G'  * the  worst  response  at  9,1)  In  the  error  region. 

9<f> 

n - the  number  of  9,0  points  used  to  define  the  feature. 
AGmax  Is  given  by: 


M - M' 


Where 


M 


Max  for  9 ,0  in  a side  lobe  or  null. 

Min  [G0(j)]  for  9,0  in  the  main  beam. 


M* 


Max 

Mln 


for  9,0  in  a side  lobe  or  null, 
for  9,0  in  the  main  beam. 


A negative  value  for  A G or  AGmax  Indicates  that  the  error  has  produced 
an  Improvement  in  the  feature.  This  will  not  occur  due  to  error  in 
and  a^,  but  may  for  error  in  1 . 

4.6  Measurement  of  Operating  Wavelength  Error  Impact 

The  designs  were  synthesized  using  an  operating  wavelength  (X)  of 
1 unit  length.  In  order  to  determine  sensitivity  to  error  in  X , 
patterns  were  made  for  each  design  where  all  element  amplitude  and  phase 
settings  were  maintained  at  the  design  solution  values,  and  operating 
wavelength  was  set  at  .98,  .99,  1.01,  and  1.02  units  of  length.  The  two 
sensitivity  measures  defined  above,  A~G  and AGmax,  were  then  measured  on 
each  feature. 
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4.7  Results 

Figures  4.12  to  4.45  are  contour  plots  of  patterns  1 through  5, 
showing  the  maximum  impacts  of  1 and  2 percent  error  In  amplitude  and 
phase  settings  on  the  features  outlined. 

Figures  4.46  to  4.50  are  graphs  of  feature  sensitivity  to  error  In 
A^  and  measured  In  terms  of  A<5  for  each  pattern.  Figures  4.51  to 
4.55  are  graphs  of  feature  sensitivity  to  error  in  A^  and  ct^,  measured 
in  terms  of  AGmax  for  each  pattern. 

Figures  4.56  and  4.57  are  graphs  of  main  beam  sensitivity  to  error 
In  A^  and  a^,  across  the  5 patterns,  measured  In  terms  of  AG,  and  AGmax 
respectively.  Figures  4.58  and  4.59  are  graphs  of  aggregated  side  lobe 
sensitivity  to  error  In  A^  and  across  the  5 patterns,  measured  in 
terms  of  A3,  and  AGmax  respectively. 

Tables  4«S  to  4.39  show  the  amplitude  and  phase  values  for  each 
pattern  generated  by  MAXIM. 

Table  4.40  summarizes  the  data  for  amplitude  and  phase  error 
sensitivity.  Table  4«42  presents  the  data  for  operating  wavelength 
sensitivity.  As  the  figures  indicate  that  sensitivity  to  £ 2 percent 
absolute  error  in  \ Is  quite  low,  it  was  felt  that  little  Information 
was  to  be  gained  by  the  presentation  of  the  associated  contour  plots. 

C 


H 


.. 


I 


49 


1 


aoaisiost  »t  10,10 
ran*  au 


0 

5 10 

5 

70  25 

10  15  10 

5 

0 55  60  65 

10 

15  80  95  90  95  00  05 

0 

5 

20  25 

10 

IS  40  45  50  55 

0 

15 

10 

15 

10 

0 

•31131131 

31 

131 131131131121131 

31131131131131 

131131131131111131131131131131131131131131131131131111131131131111 

5 

58 

58 

58 

59  56 

59 

58 

56 

58 

57 

57 

57 

57  57 

57 

57 

57 

57 

57 

57 

57 

57 

48 

59 

56 

58 

59 

56 

58 

58 

58 

10 

52 

52 

51 

51  20 

5C 

4« 

49 

48 

47 

46 

45  45 

45 

45 

45 

45 

46 

46 

47 

49 

49 

49 

49 

50 

53 

51 

51 

52 

52 

15 

41 

61 

58 

56  54 

52 

48 

46 

45 

43 

42 

39  39 

39 

36 

39 

39 

41 

42 

43 

45 

46 

49 

50 

52 

54 

56 

59 

60 

• 1 

20 

50 

50 

50 

51  51 

• | 

57 

54 

49 

45 

42 

36  35 

34 

34 

35 

36 

39 

42 

45 

49 

54 

57 

55 

53 

51 

51 

so 

50 

50 

25 

40 

40 

41 

42  43 

43 

44 

44 

46 

49 

46 

35  33 

32 

32 

33 

35 

41 

46 

48 

46 

44 

44 

43 

43 

43 

•2 

41 

40 

40 

30 

11 

31 

32 

33  15 

37 

39 

38 

37 

37 

36 

35  33 

31 

31 

33 

35 

39 

38 

37 

37 

36 

39 

39 

37 

35 

33 

32 

31 

31 

15 

28 

28 

29 

29  30 

32 

36 

36 

34 

32 

32 

32  31 

30 

30 

31 

32 

32 

32 

32 

34 

36 

36 

34 

32 

30 

29 

20 

26 

29 

40 

32 

32 

31 

31  21 

32 

17 

36 

34 

30 

29 

28  28 

27 

27 

28 

26 

28 

29 

30 

34 

36 

37 

34 

32 

31 

31 

31 

32 

32 

15 

35 

35 

34 

34  33 

34 

43 

45 

35 

30 

28 

27  27 

27 

27 

27 

27 

27 

28 

30 

35 

45 

43 

36 

34 

33 

3« 

34 

35 

It 

50 

33 

33 

31 

31  30 

30 

30 

n 

12 

32 

30 

29 

31  32 

34 

34 

32 

31 

30 

29 

30 

32 

32 

31 

30 

30 

30 

31 

31 

32 

33 

55 

64 

63 

46 

40  36 

34 

31 

30 

29 

29 

31 

30  29 

31 

31 

29 

30 

37 

11 

29 

24 

30 

31 

32 

34 

36 

40 

46 

53 

64 

50 

31 

31 

32 

32  34 

36 

47 

35 

31 

30 

3) 

27  30 

35 

35 

30 

27 

44 

33 

10 

31 

35 

47 

44 

36 

34 

32 

32 

31 

31 

45 

10 

29 

28 

28  29 

29 

32 

37 

36 

34 

33  26 

18 

18 

26 

33 

32 

34 

33 

36 

37 

32 

30 

29 

29 

20 

29 

29 

29 

38 

70 

13 

33 

31 

30  30 

30 

32 

33 

35 

30 

26  15 

10 

10 

15 

26 

30 

27 

30 

33 

35 

33 

32 

32 

30 

30 

30 

31 

32 

33 

75 

42 

41 

37 

36  36 

40 

36 

35 

36 

28 

19  10 

• 

" 

6 

10 

19 

27 

28 

33 

36 

35 

36 

65 

40 

38 

36 

37 

39 

42 

40 

68 

66 

56 

50  43 

37 

31 

31 

37 

30 

11  7 

• 

3 

7 

13 

27 

30 

41 

37 

31 

31 

33 

37 

43 

50 

56 

62 

68 

05 

60 

54 

40 

36  33 

31 

29 

29 

29 

31 

36 

T 

2 

2 

4 

9 

25 

36 

34 

31 

29 

29 

29 

31 

33 

36 

40 

46 

60 

90 

126 

56 

37 

31  30 

28 

29 

29 

29 

41 

7U 

1 

1 

3 

7 

23 

41 

31 

29 

29 

29 

28 

28 

30 

33 

37 

44 

56126 

95 

60 

54 

40 

36  33 

31 

29 

29 

31 

36 

9 4 

2 

1 

2 

4 

9 

25 

36 

34 

31 

29 

29 

29 

31 

33 

36 

40 

46 

54 

60 

100 

68 

64 

56 

50  43 

37 

21 

31 

37 

30 

13  7 

3 

3 

7 

13 

27 

30 

41 

37 

31 

31 

33 

37 

43 

50 

56 

62 

66 

60 

105 

42 

41 

37 

36  34 

40 

38 

35 

36 

28 

19  10 

6 

6 

10 

19 

27 

28 

33 

36 

35 

39 

65 

40 

36 

36 

37 

39 

41 

42 

no 

13 

33 

31 

30  30 

30 

32 

3 J 

35 

30 

26  15 

10 

10 

15 

26 

27 

10 

33 

35 

33 

32 

32 

30 

3C 

30 

31 

32 

33 

33 

115 

)0 

29 

29 

29  28 

29 

32 

37 

36 

34 

33  26 

18 

18 

26 

33 

12 

34 

13 

36 

37 

32 

30 

29 

28 

28 

28 

29 

29 

30 

120 

31 

31 

32 

32  34 

36 

47 

35 

31 

33 

27  30 

35 

35 

30 

27 

44 

33 

30 

31 

35 

47 

44 

36 

34 

32 

32 

31 

31 

31 

125 

64 

61 

46 

40  36 

34 

31 

30 

29 

29 

31 

30  29 

31 

31 

29 

30 

37 

31 

29 

29 

30 

31 

32 

34 

36 

40 

46 

53 

63 

64 

130 

33 

33 

31 

31  30 

3C 

31 

32 

32 

24 

31  32 

34 

34 

32 

31 

30 

29 

30 

32 

32 

31 

30 

30 

30 

31 

31 

32 

33 

33 

135 

15 

35 

34 

34  33 

34 

43 

45 

35 

28 

27 

27  27 

27 

27 

27 

27 

27 

28 

30 

35 

45 

43 

36 

34 

33 

34 

34 

35 

35 

35 

1«0 

12 

32 

31 

31  31 

32 

17 

38 

34 

29 

28  28 

27 

27 

26 

28 

28 

29 

30 

34 

38 

37 

14 

32 

31 

31 

31 

32 

32 

32 

145 

28 

28 

29 

29  30 

32 

36 

36 

34 

32 

32  31 

30 

30 

11 

32 

32 

32 

32 

34 

36 

36 

34 

32 

30 

29 

28 

29 

29 

28 

150 

11 

31 

32 

33  35 

37 

39 

39 

37 

39 

35  33 

31 

31 

33 

35 

39 

38 

17 

37 

36 

39 

39 

37 

35 

33 

32 

31 

31 

31 

155 

40 

40 

40 

41 

42  43 

43 

41 

44 

46 

46 

35  33 

32 

32 

13 

35 

41 

46 

48 

46 

44 

43 

43 

43 

43 

42 

41 

40 

40 

48 

140 

50 

5C 

50 

51  51 

53 

57 

54 

■9 

42 

36  35 

34 

34 

35 

36 

39 

42 

45 

49 

54 

57 

55 

S3 

51 

51 

50 

50 

50 

50 

145 

61 

61 

60 

58 

56  54 

52 

sc 

4i 

46 

45 

42 

30  39 

36 

38 

39 

39 

40 

41 

42 

4 J 

45 

46 

48 

50 

52 

54 

56 

58 

60 

61 

61 

170 

52 

52 

51 

51  50 

50 

49 

49 

49 

46 

46 

45  45 

•5 

45 

45 

45 

46 

46 

47 

48 

48 

49 

49 

50 

50 

51 

51 

52 

52 

52 

•75 

58 

58 

58 

59  56 

50 

58 

58 

59 

57 

57 

57  57 

57 

57 

57 

57 

57 

57 

57 

57 

58 

59 

58 

58 

58 

58 

59 

50 

SO 

58 

180 

131131131 

111 

131131131131131131 

i»za  ctottciTXca  zi 

•8.146 

31131131111111 

1311311311311311 31 1311311111311 311311 31131 131131 1311 31 1311 311 311 11 

f 

I 

1 

J 


1 


r 


i 

II 


i' 

1 


2 


Figure  4.12 
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Pattern  l.  1 Percent  Error 
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Figure  4.18 

Pattern  1.  1 Percent  Errot 

Maximum  Impact  on  the  Null 
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Figure  4.21 

Pattern  2.  2 Percent  Error 

Maximum  Impact  on  the  Main  Beam 
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Pattern  2.  2 Percent  Error 

Maximum  Impact  on  Side  Lobe  2 
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Figure  4.26 

Pattern  3.  1 Percent  Error 

Maximum  Impact  on  the  Main  Beam 
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Figure  4.32 

Pattern  3.  1 Percent  Error 

Maximum  Impact  on  the  Null 
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Figure  A. 36 

Pattern  A.  1 Percent  Error 
Maximum  Impact  on  Side  Lobe  1 
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Figure  4.45 

Pattern  5.  2 Percent  Error 

Maximum  Impact  on  Side  Lobe  2 
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Table  4 .6 


Pattern  1 Amplitude  and  Phase  Settings 


Maximum  Impact  on  the  Main  Beam.  1 percent  error. 


Element  Number 

Ai 

“i 

1 

0.99001 

41.49132 

2 

0.81989 

46.88094 

3 

0.54476 

57.98431 

4 

0.29153 

77.83815 

5 

0.13375 

121.66500 

6 

0.74183 

25.63405 

7 

0.56344 

37.32784 

8 

0.37299 

43.03300 

9 

0.18821 

82.60611 

10 

0.01031 

-19.88239 

11 

0.48664 

8.66110 

12 

0.32213 

7.94737 

13 

0.16964 

10.60483 

14 

0.05548 

-69.35330 

15 

0.02564 

-107.34150 

16 

0.36302 

-23.98924 

17 

0.23218 

-54.83496 

18 

0.17712 

-54.10857 

19 

-0.00683 

-100.55660 

20 

0.04657 

-91.89992 

21 

0.18199 

-75.43584 

22 

0.15151 

-53.00188 

23 

0.08324 

-112.88140 

24 

0.01772 

-103.04570 

25 

0.00364 

-155.93350 
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Table  4.7 


Pattern  1.  Amplitude  and  Phase  settings. 


Maximum  Impact  on 

the  Main  Beam. 

2 Percent  Error. 

Element  Number 

Ai 

ai 

1 

0.98007 

37.91867 

2 

0.80996 

43.30829 

3 

0.53490 

54.41331 

4 

0.30084 

74.26247 

5 

0.12382 

118.08930 

6 

0.75111 

22.06139 

7 

0.57337 

33.75526 

8 

0.38292 

39.46039 

9 

0.19814 

79.03043 

► 

10 

0.02023 

-16.30980 

11 

0.49657 

12.23315 

12 

0.33206 

11.51942 

13 

0.17957 

14.17688 

1 

14 

0.04555 

-65.77925 

15 

0-01571 

-103.76740 

16 

0.37295 

-20.41682 

17 

0.22225 

-51.26228 

I 

18 

0.18361 

-50.53598 

19 

-0.01676 

-103.90100 

20 

0.03664 

-88.32588 

f* 

21 

0.17206 

-71.86179 

V 

22 

0.14157 

-49.42924 

\ 

23 

0.07331 

-109.30740 

A 

■ i 

24 

0.00779 

-99.47165 

1 

25 

-0.00629 

-153.67570 

1- 
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Table  4.8 


Pattern  1.  Amplitude  and  Phase  settings. 
Maximum  Impact  on  Side  lobe  1.  1 Percent  Error. 


Element  Number 

Ai 

1 

1.00994 

48.68059 

2 

0.83987 

54.07021 

3 

0.54873 

65.17307 

4 

0.27155 

85.02690 

5 

0.15368 

128.85380 

6 

0.72184 

25.63399 

7 

0.54346 

37.32780 

8 

0.35301 

43.03294 

9 

0.18796 

82.72293 

10 

0.01031 

-19.90404 

11 

0.46666 

8.66104 

12 

0.30215 

7.94732 

13 

0.14966 

10.60477 

14 

0.07143 

-76.54237 

15 

0.04562 

-114.53050 

16 

0.36302 

-31.17850 

17 

0.25216 

-62.02150 

18 

0.19710 

-61.29488 

19 

0.01315 

-93.37694 

20 

0.06655 

-98.86388 

21 

0.20197 

-82.62491 

22 

0.17149 

-60.18817 

23 

0.08324 

-120.07050 

24 

0.01772 

-110.23470 

25 

0.01760 

-156.83910 
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Table  4.9 


Pattern  1.  Amplitude  and  Phase  settings. 
Maximum  Impact  on  Side  Lobe  1.  2 Percent  Error. 


Element  Number 


1.01985 

0.84979 

0.53490 

0.26162 

0.12382 

0.71192 

0.53353 

0.36832 

0.19814 

-0.01318 

0.45673 

0.29222 

0.13973 

0.08538 

0.05555 

0.37295 

0.26209 

0.20703 

0.02307 

0.07647 

0.21190 

0.18141 

0.07331 

0.00779 

0.02713 


52.25250 

57.63838 

68.74329 

88.59712 

132.42400 

22.06139 

33.75526 

39.46039 

89.43881 

-19.15211 

12.23315 

11.51942 

14.17688 

-80.11287 

-118.10100 

-34.75050 

-65.59200 

-64.86537 

-100.20710 

-102.61650 

-86.19540 

-63.75867 

-123.64100 

-113.80520 

-160.40960 
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Table  4 .10 

Pattern  1.  Amplitude  and  Phase  settings. 
Maximum  Impact  on  Side  Lobe  2.  1 Percent  Error. 


Element  Number 

A1 

“i 

1 

1.00961 

41.74815 

2 

0.82021 

53.95192 

3 

0.56436 

58.09955 

4 

0.29120 

84.90733 

5 

0.13407 

128.73420 

6 

0.74150 

32.70503 

7 

0.54379 

44.39886 

8 

0.37266 

43.15051 

9 

0.18788 

89.67531 

10 

-0.00059 

-26.80592 

11 

0.48631 

8.54297 

12 

0-30248 

0.87618 

13 

0.16931 

3.53352 

14 

0.05580 

-76.27386 

15 

0.04529 

-114.38990 

16 

0.36269 

-24.10683 

17 

0.24270 

-61.90031 

18 

0.19677 

-54.22615 

19 

-0.00651 

-93.50945 

20 

0.06622 

-98.96780 

21 

0.18231 

-75.55025 

22 

0.15183 

-60.06697 

23 

0.10313 

-112.99580 

24 

0.01804 

-110.11350 

25 

0.02329 

-149.76440 
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Table  4.11 


Pattern  1.  Amplitude  and  Phase  settings. 


Maximum  Impact 

on  Side  Lobe  2. 

2 Percent  Error. 

it  Number 

A1 

“i 

1 

1.01988 

52.25273 

2 

0.80996 

57.64056 

3 

0.57462 

54.41354 

4 

0.30146 

88.59930 

5 

0.12382 

132.42620 

6 

0.75175 

36.39546 

7 

0.53353 

33.75549 

8 

0.38292 

39.46065 

9 

0.19814 

93.36726 

10 

0.02024 

-30.64375 

11 

0.49657 

12.23329 

12 

0.29222 

-2.81394 

13 

0.17957 

-0.15656 

14 

0.04555 

-80.11456 

15 

0.05555 

-103.83400 

16 

0.36733 

-20.41698 

17 

0.26209 

-65.59369 

18 

0.20703 

-50.53621 

19 

-0.01676 

-93.19521 

20 

0.07647 

-100.87590  1 

21 

0.17206 

-71.86348  r- 

22 

0.18135 

-63.76036 

23 

0.07893 

-109.30900 

24 

0.00779 

-113.80690 

25 

0.03355 

-146.07760 
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Table  4.12 


Pattern  1.  Amplitude  and  Phase  settings. 
Maximum  Impact  on  the  Null.  1 Percent  Error. 


Element  Number 

A1 

a 

i 

1 

0.99021 

41.49115 

2 

0.81994 

52.09497 

3 

0.56463 

65.17279 

4 

0.29153 

85.02663 

5 

0.15373 

128.85350 

6 

0.72342 

32.82320 

7 

0.54346 

44.51698 

8 

0.35323 

43.03284 

9 

0.16823 

82.72266 

10 

0.01031 

-27-07143 

11 

0.48664 

1.47198 

12 

0.30764 

0.75828 

13 

0.14966 

10.60474 

14 

0.05548 

-76.54156 

15 

0.02570 

-114.52970 

16 

0.34304 

-23-98907 

17 

0.23218 

-62.02069 

18 

0.19710 

-59.32150 

19 

0.01315 

-93.49113 

20 

0.06655 

-92.01680 

21 

0.20197 

-75.55272 

22 

0.17149 

-60.11172 

23 

0.08329 

-120.06970 

24 

0.01772 

-110.23390 

25 

0.02362 

-156.83830 
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Table  4.13 

Pattern  1.  Amplitude  and  Phase  settings. 
Maximum  Impact  on  the  Mull  . 2 Percent  Error. 


Element  Number 

A1 

ai 

1 

0.98647 

37.91867 

2 

0.81008 

53.71478 

3 

0.57450 

68.74329 

4 

0.30146 

88.59712 

5 

0.16365 

132.42400 

6 

0.71232 

36.39523 

7 

0.53353 

48.08904 

8 

0.34950 

39.46039 

9 

0.15830 

79.03043 

10 

0.02023 

-30.64362 

11 

0.49657 

-2.10023 

12 

0.30314 

-2.81393 

13 

0.13973 

14.17688 

14 

0.04555 

-80.11287 

15 

0.01583 

-118.10100 

16 

0.33311 

-20.41682 

17 

0.22225 

-65.59200 

18 

0.20703 

-60.94119 

19 

0.02307 

-89.80022 

20 

0.07647 

-88.32588 

21 

0.21190 

-71.86179 

22 

0.18141 

-61.45319 

23 

0.07343 

-123.64100 

24 

0.00779 

-113.80520 

25 

0.03355 

-160.40960 
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Table  4 .14 

Pattern  2.  Amplitude  and  Phase  settings. 
Maximum  Impact  on  the  Main  Beam.  1 Percent  Error. 


Element  Number 

Ai 

°i 

1 

0.99001 

82.99451 

2 

0.71493 

72.69763 

3 

0.41100 

42.82335 

4 

0.22526 

-13.39055 

5 

0.15096 

-55.64073 

6 

0.66943 

78.32292 

7 

0.50047 

59.45757 

8 

0.24733 

17.08455 

9 

0.25547 

-29.58571 

10 

0.12430 

-48.79063 

11 

0.25718 

55-11913 

12 

0.22163 

16.72641 

13 

0.28098 

-5.35801 

14 

0.22563 

-32.30592 

15 

0.12461 

-58.47069 

16 

0.15506 

-27.54604 

17 

0.15947 

-25.03529 

18 

0.22672 

-33.80103 

19 

0.13454 

-50.83426 

20 

0-01818 

-64.63976 

21 

0.12854 

-62.08603 

22 

0.10943 

-50.34950 

23 

0.05815 

-44.35547 

24 

0.06849 

-60.85086 

25 

-0.00923 

-81.58513 
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Table  4.15 


Pattern  2*  Amplitude  and  Phase  settings. 
Maximum  Impact  on  the  Main  Beam.  2 Percent  Error. 


Element  Number 

Ai 

°i 

1 

0.98008 

79.42007 

2 

0.70500 

69.12321 

3 

0.42092 

39.25085 

4 

0.23519 

-9.81820 

5 

0.14103 

-52.06812 

6 

0.65950 

74.74850 

7 

0.51040 

55.88519 

8 

0.25725 

15.54277 

9 

0.26540 

-26.01320 

10 

0.12097 

-45.21815 

11 

0.26710 

51.54634 

12 

0.23156 

14.19118 

13 

0.29091 

-1.78578 

14 

0.23555 

-28.73344 

15 

0.11468 

-54.89915 

16 

0.16498 

-23.97356 

17 

0.16940 

-21.46283 

18 

0.23665 

-30.22852 

19 

0.12461 

-47.26175 

20 

0.00825 

-61.06621 

21 

0.11861 

-58.51248 

22 

0.11317 

-46.77702 

23 

0.06899 

-40.78296 

24 

0.05857 

-57.27869 

25 

-0.01916 

-85.18558 
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Table  4.16 


Pattern  2.  Amplitude  and  Phase  settings. 
Maximum  Impact  on  Side  Lobe  1.  1 percent  Error. 


Element  Number 

Ai 

°i 

1 

0.99001 

82.99374 

2 

0.71493 

72.69687 

3 

0.39101 

42.82324 

4 

0.22526 

-20.57953 

5 

0.17094 

-62.80807 

6 

0.68941 

85.51195 

7 

0.50047 

66.64658 

8 

0.22735 

22.23720 

9 

0.23549 

-29.58563 

10 

0.12430 

-55.97975 

11 

0.23719 

55.11902 

12 

0.22163 

15.79092 

13 

0.28098 

-12.54722 

14 

0.22563 

-39.49507 

15 

0.14459 

-62.30493 

16 

0.13508 

-27.54596 

17 

0.13949 

-25.03519 

18 

0.20674 

-33.80095 

19 

0.13454 

-58.02119 

20 

0.01818 

-71.82802 

21 

0.14852 

-62.08505 

22 

0.12908 

-57.53702 

23 

0.07780 

-51.54457 

24 

0.08847 

-61.31598 

25 

0.01075 

-74.48824 

Pattern  2.  Amplitude  and  Phase  settings. 
Maximum  Impact  on  Side  Lobe  1.  2 Percent  Error. 


Element  Number 

A1 

ai 

1 

0.98007 

79.41789 

2 

0.70500 

69.12102 

3 

0.38108 

39.25060 

4 

0.23518 

-24.15146 

5 

0.18087 

-66.35463 

6 

0.69933 

89.08096 

7 

0.51039 

70.21561 

8 

0.21741 

25.80904 

9 

0.22556 

-26.01297 

10 

0.11437 

-59.54799 

11 

0.22726 

51.54608 

12 

0.23156 

12.21849 

13 

0.29091 

-16.11913 

14 

0.23555 

-43.06696 

15 

0.15451 

-65.30235 

16 

0.12515 

-23.97336 

17 

0.12956 

-21.46260 

18 

0.19681 

-30.22830 

19 

0.12461 

-61.59137 

20 

0.00825 

-75.39813 

21 

0.15845 

-66.00323 

22 

0.13901 

-61.10658 

23 

0.08773 

-55.11650 

24 

0.09840 

-67.68253 

25 

0.02067 

-71.09834 
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Table  4.18 

Pattern  2. 

Amplitude  and  Phase 

settings. 

Maximum  Impact 

on  Side  Lobe  2. 

1 

Percent  Error. 

Element  Number 

Ai 

“i 

1 

1.00995 

83.11012 

2 

0.71493 

79.88518 

3 

0.41099 

42.82304 

4 

0.20528 

-13.39043 

5 

0.17094 

-62.82626 

6 

0.66943 

85.51047 

7 

0.50047 

59.57318 

8 

0.22734 

15.04794 

9 

0.25547 

-36.77461 

I 10 

0.13350 

-48.79034 

11 

0.25717 

55.19592 

12 

0.20165 

15.79082 

13 

0.28098 

-5.35789 

i 14 

0.20837 

-32.30568 

15 

0.14459 

-65.65773 

16 

0.15505 

-34.65994 

i 17 

0.13949 

-32.22432 

18 

0.22672 

-33.80075 

1 19 

0.13454 

-58.02000 

20 

0.01818 

-71.82681 

21 

0.14852 

-69.27309 

22 

0.12908 

-50.34926 

23 

0.05782 

-51.54437 

24 

0.08847 

-68.03792 

, 25 

-0.00923 

-81.67628 

( 
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Table  4.19 


Pattern  2.  Amplitude  and  Phase  settings. 
Maximum  Impact  on  Side  Lobe  2.  2 Percent  Error. 


Element  Number 

Ai 

ai 

1 

1.01983 

79.41620 

2 

0.70500 

83.45398 

3 

0.42092 

39.25041 

4 

0.19535 

-9.81802 

5 

0.18086 

-66.39648 

6 

0.65949 

89.07927 

7 

0.51039 

55.88481 

8 

0.21741 

11.47547 

9 

0.26539 

-40.34660 

10 

0.14813 

-45.21776 

11 

0.26710 

52.10831 

12 

0.19172 

12.21839 

13 

0.29090 

-1.78575 

14 

0.22138 

-28.73311 

15 

0.15451 

-69.22795 

16 

0.16498 

-37.74136 

17 

0.12956 

-35.79630 

18 

0.23664 

-30.22812 

19 

0.12461 

-61.59027 

20 

0.00825 

-75.39703 

21 

0.15845 

-72.84331 

22 

0.13901 

-46.77669 

23 

0.04789 

-55.11636 

24 

0.09840 

-71.60814 

25 

-0.01916 

-85.24651 

f 

i 

97 

A 

f 
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( 

Table  4.20 

Pattern  3. 

Amplitude  and  Phase 

settings. 

Maximum  Impact 

on  the  Main  Beam.  1 

Percent  Error. 

Clement  Number 

Ai 

“i 

1 

0.99033 

83.10712 

2 

0.71665 

72.81026 

3 

0.35546 

47.94011 

4 

0.24539 

-18.50732 

5 

0.15889 

-40.51294 

6 

0.68453 

78.43555 

7 

0.48179 

59.57018 

8 

0.25729 

16.45805 

9 

0.24596 

-27.14305 

10 

0.14811 

-57.90044 

11 

0.25928 

55.23537 

12 

0.20170 

17.95587 

13 

0.27051 

-5.47562 

14 

0.21723 

-37.42282 

15 

0.09138 

-48.98233 

16 

0.15109 

-31.51103 

i 

17 

0.15355 

-22.71280 

, 

18 

0.21828 

-33.91832 

19 

0.12946 

-45.95184 

i 

20 

0.04597 

-54.51227 

21 

0.10977 

-65.39644 

: 

22 

0.10497 

-53.82712 

23 

0.07493 

-40.11311 

24 

0.06588 

-47.36270 

4 

> 

25 

( 

-0.00894 

-76.95766 

1 

i ^ 

i- 

I 

t 

\ 

* 

i 

I i 

1 

7^3? -?j  ■ > 

* * * .JIUJL.  L V1.',. 

98 


Table  4.21 


Pattern  3.  Amplitude  and  Phase  settings. 
Maximum  Impact  on  the  Main  Beam.  2 Percent  Error. 


Element  Number 

Ai 

ai 

1 

0.98008 

79.41866 

2 

0.70639 

69.12178 

3 

0-36572 

44.25023 

4 

0.25565 

-14.81753 

5 

0.16914 

-36.82297 

6 

0.67427 

74.74707 

7 

0-49204 

55.88377 

8 

0.26755 

15.53319 

9 

0.25621 

-23.45320 

10 

0.13786 

-54.21576 

11 

0.26953 

51.54518 

12 

0.21196 

15.02099 

13 

0.28077 

-1.78571 

14 

0.22749 

-33.73288 

15 

0.08113 

-45.29236 

16 

0.16135 

-27.82106 

17 

0.16381 

-19.02280 

18 

0.22854 

-30.22832 

19 

0.11921 

-42.26193 

20 

0.03571 

-50.82227 

21 

0.09952 

-61.70944 

22 

0.09588 

-50.13727 

23 

0.07958 

-36.42317 

24 

0.05563 

-43.67273 

25 

-0.01919 

-80.56252 
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Table  4.22 

Pattern  3*  Amplitude  and  Phase  settings. 
Maximum  Impact  on  Side  Lobe  1.  1 Percent  Error. 


Element  Number 

A1 

“i 

1 

0.99001 

90.18063 

2 

0.71632 

79.88376 

3 

0.33603 

55.01173 

4 

0.22574 

-18.38959 

5 

0.13924 

-40.39520 

6 

0.70418 

85.48949 

7 

0.48212 

66.64369 

8 

0.23764 

15.04776 

9 

0.24629 

-34.21460 

10 

0.16777 

-64.97406 

11 

0.23963 

55.11768 

12 

0.18205 

17.83820 

13 

0.25108 

-5.35782 

14 

0.19758 

-37.30504 

15 

0.09105 

-48.86459 

16 

0.13144 

-31.39325 

17 

0.13390 

-22.59503 

18 

0.21861 

-40.98973 

19 

0.14912 

-53.02333 

20 

0.06562 

-61.58011 

21 

0.12943 

-72.47006 

22 

0.12463 

-60.89513 

23 

0.05528 

-39.99532 

24 

0.06556 

-47.24496 

25 

-0.00926 

-77.05322 

1 

f 

( 

Table  4.23 

100 

Pattern  3. 

Amplitude  and  Phase 

settings. 

Maximum  Impact 

on  Side  Lobe  1.  2 

Percent  Error. 

Element  Number 

A1 

“i 

1 

0.98007 

93.75113 

2 

0.70639 

83.45425 

3 

0.33230 

58.57985 

4 

0.21581 

-14.81742 

5 

0.12931 

-36.82277 

6 

0.71411 

89.03709 

7 

0.49204 

70.21419 

8 

0.22771 

11.47540 

9 

0.25621 

-37.78667 

10 

0.17770 

-68.54538 

11 

0.22970 

51.54495 

12 

0.17212 

14.26602 

13 

0.24735 

-1.78570 

14 

0.18765 

-33.73273 

15 

0.08113 

-45.29216 

16 

0.12151 

-27.82083 

17 

0.12397 

-19.02257 

18 

0.22854 

-44.56186 

19 

0.15904 

-56.59546 

20 

0.07554 

-65.15143 

21 

0.13935 

-76.04137 

* 

22 

0.13455 

-64.46645 

23 

0.04536 

-42.61943 

24 

0.05563 

-43.67253 

25 

C 

I 

-0.01919 

-79.54712 

1 i 
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Table  4.24 


Pattern  3.  Amplitude  and  Phase  settings. 
Maximum  Impact  on  Side  Lobe  2*  1 Percent  Error. 


Element  Number 

A1 

ai 

1 

1.00995 

90.18057 

2 

0.71632 

72.81178 

3 

0.33581 

55.01170 

4 

0.24572 

-25.57875 

5 

0.15922 

-47.58437 

6 

0.70418 

85.50900 

7 

0.46214 

59.57170 

8 

0.23764 

22.23676 

9 

0.24629 

-34.21460 

10 

0.16777 

-63.57092 

11 

0.25961 

62.30411 

12 

0.18205 

17.83817 

13 

0.27084 

-5.35781 

14 

0.19758 

-44.49431 

15 

0.11104 

-56.05379 

16 

0.15142 

-31.39322 

17 

0.13390 

-29.78421 

18 

0.21861 

-33.80057 

19 

0.12913 

-53-02332 

20 

0.06552 

-61.57962 

21 

0.12943 

-71.98157 

22 

0.10470 

-60.89464 

23 

0.06820 

-39.99530 

24 

0.08554 

-54.43413 

25 

-0.00926 

-69.98135 
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Table  A. 25 


Pattern  3.  Amplitude  and  Phase  settings. 
Maximum  Impact  on  Side  Lobe  2.  2 Percent  Error. 


Element  Number 

Ai 

ai 

1 

1.01987 

93.75259 

2 

0.70639 

69.12108 

3 

0.32588 

58.58133 

4 

0.25565 

-29.15099 

5 

0.16914 

-51.03870 

6 

0.71411 

89.08102 

7 

0.45221 

55.88370 

8 

0.22771 

25.80901 

9 

0.25621 

-37.78683 

10 

0.17770 

-68.54456 

11 

0.26953 

65.87613 

12 

0.17213 

14.26612 

13 

0.28077 

-1.78570 

14 

0.18765 

-48.06648 

15 

0.12096 

-59.62334 

16 

0.16135 

-27.82098 

17 

0.12397 

-33.35643 

18 

0.22854 

-30.22826 

19 

0.11920 

-56.59561 

20 

0.06434 

-65.15274 

21 

0.13935 

-70.39001 

22 

0.12109 

-64.46776 

23 

0.05657 

-36.42313 

24 

0.09547 

-58.00421 

25 

-0.01919 

-66.29224 
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Table  4.26 


Pattern  3*  Amplitude  and  Phase  settings. 


Maximum  Impact  on  the  Null 


1 Percent  Error. 


Element  Number 

Ai 

ai 

1 

0.99033 

83.10838 

2 

0.71665 

79.76546 

3 

0.33614 

54.89377 

4 

0.24539 

-18.50745 

5 

0.15889 

-40.51306 

6 

0.70385 

85.39075 

7 

0.48179 

59.57143 

8 

0.25729 

15.16562 

9 

0.22663 

-34.09666 

10 

0.14812 

-64.85544 

11 

0.25928 

55.23553 

12 

0.18238 

24.90955 

13 

0.25119 

-5.47570 

14 

0.21723 

-37.42291 

15 

0.11071 

-48.98245 

16 

0.13176 

-38.46460 

17 

0.15355 

-22.71289 

18 

0.21828 

-33.91841 

19 

0.12946 

-52.90540 

20 

0.04597 

-61.46149 

21 

0.12910 

-72.35143 

22 

0.11508 

-60.90169 

23 

0.05561 

-47.04689 

24 

0.06588 

-54.31622 

25 

0.01039 

-69.98114 
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Table  4.27 


Pattern  3. 

Amplitude  and 

Phase 

settings. 

Maximum  Impact 

on  the  Null 

2 

Percent  Error. 

it  Number 

Ai 

°i 

1 

0.98007 

79.41647 

2 

0.70639 

83.45425 

3 

0.32588 

58.57985 

4 

0.25565 

-25.22350 

5 

0.16914 

-44.31465 

6 

0.71411 

89.07954 

7 

0.49204 

55.88351 

8 

0.26754 

11.47540 

9 

0.21637 

-37.78667 

10 

0-13786 

-68.54538 

11 

0.26953 

51.54495 

12 

0.17212 

28.59955 

13 

0.24093 

-1.78570 

14 

0.22749 

-33.73273 

15 

0.12096 

-45.29216 

16 

0.12151 

-42.15465 

17 

0.16381 

-24.37172 

18 

0.22854 

-30.22809 

19 

0.11920 

-56.59546 

20 

0.03571 

-65.15143 

21 

0.13935 

-76.04137 

22 

0.11554 

-62.16019 

23 

0.04536 

-50.71310 

24 

0.05563 

-58.00290 

25 

0.02064 

-76.69783 
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Table  4.28 


Pattern  4*  Amplitude  and  Phase  settings. 
Maximum  Impact  on  the  Main  Beam.  1 Percent  Error. 


Element  Number 

Ai 

°i 

1 

0.99000 

80.22981 

2 

0.57462 

71.41367 

3 

0.15240 

30.66972 

4 

0.14752 

-58.21255 

5 

0.09088 

-96.90738 

6 

0.54779 

80.69716 

7 

0.35826 

56.28606 

8 

0.24028 

9.77547 

9 

0.19048 

-21.48817 

10 

0.09229 

-73.88510 

11 

-0.00589 

-171.35440 

12 

0.15118 

2.71470 

13 

0.27176 

9.15734 

14 

0.20907 

-10.04460 

15 

0.09625 

-8.75852 

16 

0.14645 

-95.75626 

17 

0.11651 

-41.43205 

18 

0.17666 

-3.24723 

19 

0.15050 

4.94055 

20 

0.02904 

-6.78648 

21 

0.07705 

-96.84164 

22 

0.02588 

-75.27524 

23 

0.04555 

-36.44624 

24 

0.03225 

111.78390 

25 

0.02849 

85.51938 
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Table  4.29 


Pattern  4.  Amplitude  and  Phase  settings. 
Maximum  Impact  on  the  Main  Beam.  2 Percent  Error. 


Element  Number 

Ai 

“i 

1 

0.98007 

76.53784 

2 

0-56469 

67.72171 

3 

0.16233 

27.09732 

4 

0.13759 

-54.52652 

5 

0.08095 

-93.21754 

6 

0.53786 

77.00519 

7 

0.36819 

52.71344 

8 

0.25021 

13.34765 

9 

0.20040 

-17.91579 

10 

0.08236 

-70.19528 

11 

-0.01582 

-170.89520 

12 

0.16111 

6.28663 

13 

0.28168 

12.72952 

14 

0.21899 

-6.47231 

15 

0.10617 

-5.18623 

16 

0.13653 

-92.06642 

17 

0.12637 

-37.85968 

18 

0.18658 

0.32496 

19 

0.16043 

8.51276 

20 

0.03897 

-3.21438 

21 

0.06712 

-93.15181 

22 

0.01596 

-71.58540 

23 

0.05548 

-32.87387 

24 

0.04218 

108.09200 

25 

0.03842 

81.82742 
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Table  4.30 


Pattern  4.  Amplitude  and  Phase  settings. 


Maximum  Impact  on  Side  Lobe  1.  1 Percent  Error. 


Element  Number 


a 


1 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 


0.99001 

0.57462 

0.13242 

0.14752 

0.09088 

0.56778 

0.35826 

0.24029 

0.19048 

0.11227 

-0.00589 

0.15097 

0.25178 

0.18909 

0.07627 

0.14646 

0.09658 

0.17666 

0.15050 

0.02904 

0.09703 

0.04586 

0.02557 

0.03225 

0.02849 


80.23009 

78.48587 

37.85910 

-65.28502 

-103.97980 

87.76935 

63.47284 

2.56845 

-28.67743 

-80.95758 

-178.27770 

-4.47435 

9.15739 

-10.04464 

-8.75856 

-102.82870 

-41.43213 

-10.43635 

-2.24838 

-13.97578 

-96.84274 

-82.34772 

-43.63550 

118.85610 

92.59158 


Table  4.31 


Pattern  4.  Amplitude  and  Phase  settings. 
Maximum  Impact  on  Side  Lobe  1.  2 Percent  Error. 


Element  Number 


a 


i 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 


0.98007 

0.56469 

0.12249 

0.13759 

0.10301 

0.57770 

0.36819 

0.25021 

0.20040 

0.12219 

-0.01582 

0.15954 

0.24185 

0.17916 

0.06634 

0.13652 

0.08665 

0.18658 

0.16043 

0.03897 

0.10695 

0.05579 

0.01564 

0.04218 

0.03842 


76.53615 

82.05467 

41.43098 

-68.85524 

-107.55000 

91.33815 

67.04164 

-0.98577 

-32.24942 

-82.19569 

-181.36070 

-8.04643 

12.72942 

-6.47224 

-5.18616 

-106.39890 

-37.85950 

-14.00842 

-5.82033 

-17.54759 

-93.15073 

-82.00270 

-47.20738 

122.42490 

96.16039 
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Table  4.32 

Pattern  4.  Amplitude  and  Phase  settings. 
Maximum  Impact  on  Side  Lobe  2>  1 Percent  Error. 


Element  Number 

Ai 

1 

1.00995 

87.30200 

2 

0.57462 

71.41594 

3 

0.13242 

37.85910 

4 

0.16751 

-58.21365 

5 

0.11086 

-103.97980 

6 

0.56778 

87.76935 

7 

0.33828 

56.28613 

8 

0.22031 

9.77551 

9 

0.19048 

-21.48825 

10 

0.11227 

-80.95758 

11 

0.01409 

-178.35520 

12 

0.15118 

-4.47435 

13 

0.25178 

9.15739 

14 

0.20907 

-10.04464 

15 

0.07627 

-15.94782 

16 

0.16644 

-102.82870 

17 

0.11651 

-48.62131 

18 

0.15668 

-10.43635 

19 

0.15050 

4.94060 

20 

0.00906 

-13.97578 

21 

0.09703 

-103.83660 

22 

0.02588 

-82.34772 

23 

0.04555 

-43.63550 

24 

0.03225 

118.83660 

25 

0.00851 

92.59158 
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Table  4.33 

Pattern  4.  Amplitude  and  Phase  settings. 
Maximum  Impact  on  Side  Lobe  2*  2 Percent  Error. 


Element  Number 

Ai 

ai 

1 

1.01985 

90.87250 

2 

0.56469 

67.72171 

3 

0.12249 

41.43115 

4 

0.17743 

-54.52652 

5 

0.12078 

-107.55110 

6 

0.57770 

61.33984 

7 

0.32835 

52.71344 

8 

0.21038 

13.34765 

9 

0.20040 

-17.91579 

10 

0.12219 

-84.52890 

11 

0.02401 

-177.61450 

12 

0.16111 

-8.04650 

13 

0.24185 

12.72952 

14 

0.21899 

-11.57497 

15 

0.06634 

-19.51967 

16 

0.17636 

-106.40000 

17 

0.12637 

-52.19350 

18 

0.14675 

-14.00849 

19 

0.16043 

7.43505 

20 

-0.00087 

-11.92710 

21 

0.10696 

-105.17790 

22 

0.03318 

-85.91902 

23 

0.05548 

-47.20755 

24 

0.04218 

122.38420 

25 

-0.00141 

91.33383 

\ 
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TABLE  4.34 

Pattern  5. 

Amplitude  and 

Phase  settings. 

Maximum  Impact 

on  the  Main  Beam.  1 Percent  Error. 

Element  Number 

Ai 

“l 

1 

0.65060 

99.26013 

2 

0.37304 

106.42760 

3 

0.00745 

101.81600 

4 

0.08857 

-82.23674 

5 

0.08640 

-45.24315 

6 

1.00998 

38.29744 

7 

0.66055 

33.42734 

8 

0.18020 

34.65767 

9 

0.07037 

-89.73033 

10 

0.13958 

-143.54060 

11 

0.82609 

6.09347 

! 12 

0.48946 

4.94162 

I 13 

0.18212 

-16.82225 

14 

0.06276 

-164.23010 

15 

0.05580 

-161.97100 

16 

0.52866 

-71.47241 

| 17 

0.37535 

-74.10629 

i is 

0.11292 

-89.93283 

S 19 

0.01421 

59.96715 

1 20 

0.09094 

65.32414 

21 

0.26640 

-97.65575 

- 

22 

0.16686 

-103.73350 

23 

0.04862 

-91.10423 

24 

0.06307 

39.77171 

25 

0.04109 

50.08115 
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TABLE  4.35 


Pattern  5.  Amplitude  and  Phase  settings. 
Maximum  Impact  on  the  Main  Beam.  2 Percent  Error. 


Element  Number 

A1 

ai 

1 

0.64067 

95.68428 

2 

0.36311 

102.85170 

3 

-0.00248 

103.27540 

4 

0.07864 

-78.66121 

5 

0.09632 

-41.67044 

6 

1.01987 

34.72470 

7 

0.67048 

29.85460 

8 

0.19013 

31.08493 

9 

0.06044 

-86.15480 

10 

0.12964 

-139.96510 

11 

0.83602 

9.66548 

12 

0.49939 

8.51364 

13 

0.19205 

-13.24991 

14 

0.05283 

-160.64750 

15 

0.04587 

-158.39550 

16 

0.51873 

-67.89690 

17 

0.34742 

-70.53078 

18 

0.10299 

-86.35732 

19 

0.02414 

56.39423 

20 

0.10087 

61.74829 

21 

0.25647 

-94.08023 

22 

0.15693 

-100.15800 

23 

0.03870 

-87.52872 

24 

0.07300 

36.19897 

25 

0.05101 

46.50844 
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TABLE  4.36 


Pattern  5.  Amplitude  and  Phase  settings. 
Maximum  Impact  on  Side  Lobe  1.  1 Percent  Error. 


Element  Number 

A1 

ai 

1 

0.67058 

99.37532 

2 

0.39302 

106.54280 

3 

0.02743 

108.96940 

4 

0.10855 

-89.40341 

5 

0.08640 

-52.43205 

6 

0.99001 

38.29704 

7 

0.64057 

33.46292 

8 

0.16022 

41.81230 

9 

0.07037 

-96.91640 

10 

0.15956 

-150.72670 

11 

0.82609 

6.09326 

12 

0.48946 

4.94142 

13 

0.18212 

-16.82198 

14 

0.06276 

-164.31280 

15 

0-07572 

-162.08570 

16 

0.54864 

-78.65848 

17 

0.37733 

-81.29237 

18 

0.13289 

-97.11890 

19 

-0.00577 

59.99971 

20 

0.09094 

72.51126 

21 

0.26640 

-97.77045 

22 

0.16685 

-103.84820 

23 

0.04862 

-91.21893 

24 

0.06307 

39.77130 

25 

0.02111 

50.08075 
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TABLE  4.37 


Pattern  5.  Amplitude  and  Phase  settings. 


Maximum  Impact  on  Side  Lobe  1.  2 Percent  Error. 


( 


Element  Number 

Ai 

“i 

1 

0.68051 

106.09150 

2 

0.40294 

113.25900 

3 

0-03736 

110.26660 

4 

0.11847 

-92.95166 

5 

0.07910 

-56.00414 

6 

0.98007 

34.72459 

7 

0.63064 

29.85451 

8 

0.15029 

43.11046 

9 

0.07766 

-100.48820 

10 

0.16948 

-154.29850 

11 

0.83602 

9.66543 

12 

0.49939 

8.51360 

13 

0-17931 

-13.24987 

14 

0.07365 

-160.62240 

15 

0.08558 

-158.39530 

16 

0.55856 

-82.23029 

17 

0.38725 

-84.86417 

18 

0.14282 

-100.69070 

19 

-0.01570 

58.69946 

20 

0.10086 

76.08180 

21 

0.25647 

-94.08002 

22 

0.15692 

-100.15770 

23 

0.03870 

-87.52849 

24 

0.07300 

44.78825 

25 

0.01378 

46.98734 
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TABLE  4.38 

Pattern  5.  Amplitude  and  Phase  settings. 
Maximum  Impact  on  Side  Lobe  2.  1 Percent  Error. 


Element  Number 

A1 

1 

0.65060 

106.44720 

2 

0.39302 

113.61470 

3 

0.00745 

101.93120 

4 

0.10855 

-89.42329 

5 

0.06642 

-45.24286 

6 

1.00995 

45.48639 

7 

0.64057 

33.42696 

8 

0.16022 

41.84660 

9 

0.07037 

-89.84552 

10 

0.15956 

-143.65580 

11 

0.80611 

-1.09571 

12 

0.48946 

4.94144 

13 

0.16214 

-24.01114 

14 

0.06276 

-171.38470 

15 

0.07578 

-165.52790 

16 

0.54864 

-78.65897 

17 

0.35735 

-81.29286 

18 

0.13290 

-90.04802 

19 

-0.00577 

60.04353 

20 

0.09094 

72.43388 

21 

0.26640 

-104.84230 

22 

0.18683 

-110.92000 

23 

0.04862 

-91.21942 

24 

0.04309 

46.96063 

25 

0.03351 

50.08078 
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TABLE  4.39 

Pattern  5.  Amplitude  and  Phase  settings. 


Maximum  Impact  on  Side  Lobe  2.  2 Percent  Error. 


Element  Number 

Ai 

“i 

1 

0.64067 

110.01770 

2 

0.40294 

117.18520 

3 

-0.00248 

104.43720 

4 

0.11847 

-92.99461 

5 

0.05649 

-41.67038 

6 

1.00713 

49.05843 

7 

0.63064 

29.85451 

8 

0.15029 

45.41866 

9 

0.06044 

-86.15459 

10 

0.16948 

-139.96490 

11 

0.79618 

-4.66772 

12 

0.49939 

8.51360 

13 

0.15221 

-27.58324 

14 

0.05283 

-174.95600 

15 

0.08570 

-162.32360 

16 

0.55856 

-82.23029 

17 

0.34742 

-84.86417 

18 

0.14282 

-86.35709 

19 

-0.01570 

58.69946 

20 

0.10086 

73.77478 

21 

0.25647 

-108.41360 

22 

0.19676 

-114.49130 

23 

0.03870 

-87.52849 

24 

0.03317 

50.53267 

25 

0.01147 

46.50833 
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Table  4.40.  Summary  of  A , a Impact  Data. 


Average 
Impact 
in  dB. 


Error  Range  In  Z 


Figure  4.49 
Pattern  4 

AG  vs • Error  Range  by  Feature. 


•,  -a . 


- 


Error  Range  In  * 


Figure  4.51 
Pattern  1 

A Gmax  vs.  Error  Range  by  Feature. 


Error  Range  in  Z 


Figure  4.53 
Pattern  3 

A Gmax  vs,  Error  Range  by  Feature. 


Error  Range  in  % 


Figure  4.54 
Pattern  4 

4 Gaax  va.  Error  Range  by  Feature. 


Error  Range  in  Z 


Figure  4.56 

Main  Bean  Impact 

AG  vs.  Error  Range  by  Pattern. 
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Error  Range  In  Z 


Figure  4*58 

Aggregated  Side  Lobe  Impact. 
AG  va.  Error  Range  by  Pattern. 
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Figure  4.59 

Aggregated  Side  Lobe  Impact. 

A Gmax  vs.  Error  Range  by  Pattern. 


Table  4.42.  Summary  of  X Error  Impact  Data. 


CHAPTER  V 


DISCUSSION  OF  THE  RESULTS 
AND  EXTENSIONS  OF  THE  METHOD 


5.1  Implications  of  the  results 

The  analysis  presented  is  meant  to  illustrate  the  power  of  the 
MAXIM-IOS  method,  and  not  to  generate  sensitivity  data  for  a particular 
design  problem.  None  the  less,  some  trends  which  may  be  tentatively 
identified  in  the  data  deserve  comment. 

1.  The  nulls  appear  as  the  most  sensitive  feature  in  the 
pattern  by  a factor  of  about  three  over  the  side  lobes. 

2.  Among  the  features  examined,  the  main  beams  showed  the 
least  sensitivity  overall.  However,  among  the  main  beams,  the 
square  or  rectangular  shapes  were  much  more  sensitive  at  the  2 
percent  error  range  than  the  round  beam. 

3«  In  general,  the  broader  the  main  beam,  the  greater  the 
side  lobe  sensitivity  appeared  to  be. 

4.  The  minimum  elevation  points  on  the  error  surface  almost 
universally  occurred  at  the  edge  of  the  error  Interval  of  each 
decision  variable. 

The  first  three  observations  accord  well  with  intuitive 
expectations  for  feature  sensitivity.  The  fourth  observation  is 
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I exactly  what  would  be  expected  for  a well  behaved  error  surface.  That 

\ 

Is,  the  error  regions  appear  to  be  tightly  enough  constrained  about  an 
optimum  to  contain  only  one  local  minimum*  It  Is  to  be  expected  that 
such  a minimum  will  lie  on  an  edge,  as  far  from  the  maximum  as  possible. 
Thus  It  can  be  said  that  for  a 10  by  10  array  at  least,  MAXIM-generated 
sensitivity  values  can  be  accepted  with  some  confidence  within  the  two 
percent  error  region. 


( 


5.2  Potential  uses  for  MAXIM 

The  worst  performance  possslble  for  an  array  design  or  feature, 
deriving  from  error  Is,  in  itself,  a valuable  piece  of  information. 
However,  MAXIM  can  be  used  to  generate  other  types  of  Information  as 
well.  In  Chapter  four,  only  one  type  of  experiment  which  may  be 
designed  around  MAXIM  was  described.  There,  a method  to  determine  the 
relationship  between  directivity  pattern  and  sensitivity  was  examined. 
If  such  an  experiment  were  to  be  performed  on  a series  of  alternative 
pattern  designs,  the  procedure  could  be  refined  somewhat.  Presumably, 
the  alternatives  to  be  evaluated  in  an  actual  application  are  not  likely 
to  differ  from  each  other  aa  radically  as  the  sample  in  the  above  study. 
The  sensitivity  of  a pattern  might,  for  instance,  be  examined  as  one 
feature  in  a pattern  changes  over  some  Interval.  It  may  prove  that  In 
some  cases.  It  is  possible  to  greatly  decrease  sensitivity  through  small 
modifications  in  the  design  specifications.  It  should  be  noticed  that 
in  such  a case  the  IOS  will  tend  to  be  based  on  a meaningful  set  of 
standards. 

It  would  also  be  of  value  to  take  MAXIM  values  on  a finer  set  of 
error  ranges  than  those  used  In  this  study.  The  relative  rate  of  change 
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in  maximum  degradation  of  given  features  over  an  error  range  is  of 
importance,  since  it  can  give  an  indication  of  the  expected  degradation. 
For  Instance,  suppose  two  features  were  to  be  compared  for  sensitivity, 
and  the  following  MAXIM  data  were  obtained.  At  one  percent  error  the 
first  feature  had  a sensitivity  value  of  5 dB , while  the  second  had  a 
sensitivity  value  of  10  dB.  At  two  percent  error  the  sensitivity  value 
of  both  features  was  15  dB.  In  this  case,  the  expected  performance  of 
the  first  feature  would  be  worse  than  the  second. 

Element  amplitude  and  phase  coefficient  error  ranges  could  be 
analyzed  separately,  allowing  these  ranges  to  take  on  realistic  values 
for  the  particular  system  being  studied. 

If,  in  a given  application,  array  pattern  sensitivity  is  to  be 
taken  into  account  at  the  design  stage,  it  is  possible,  by  holding  any 
two  of  the  following  three  array  characteristics  constant,  to  perform  a 
sensitivity  analysis  similar  to  the  one  in  this  study  on  the  third. 

1.  Array  size  and  geometry. 

2.  Directivity  pattern  design. 

3.  Error  ranges. 

Number  three  assumes  the  existence  of  a reliable  means  of  converting  the 
effects  of  error  sources,  such  as  those  mentioned  in  Chapter  two,  into 
error  ranges  in  A^  and  and  conversely.  Then,  sensitivity  value  vs. 
error  range  graphs,  such  as  figure  4.56  could  be  used  to  specify 
acceptable  tolerances  for  array  or  component  manufacture. 

Finally,  the  arrays  examined  in  this  study  have  all  been  of  the 
planar  symmetric  type.  MAXIM  may  easily  be  extended  to  the  more  general 
case  of  conformal  array  sensitivity  analysis,  since  NLGP  Conformal,  a 
non-linear  goal  programming  code  specifically  for  the  conformal  array 
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case,  has  already  been  developed*  A modification  of  this  code  to  yield 
a conformal  version  of  MAXIM  would  also  make  it  possible  to  analyze  non- 
symmetric  error  In  planar  symmetric  arrays. 


5.3  Summary 

The  MAXIM-IOS  method  of  sensltivlty/impact  analysis  has  been 
presented  and  Illustrated  with  a sample  study  on  five  rather  complex 
array  patterns.  The  resulting  measure,  maximum  Impact,  represents  the 
worst  realization  which  a particular  feature  of  a pattern  can  assume  due 
to  a given  amount  of  error.  The  change  in  the  value  of  this  realization 
as  error  in  amplitude  and  phase  settings  increases  is  also  an  Indicator 
of  the  shape  of  the  response  surfance  in  the  error  region. 

Since  MAXIM  is  based  on  an  optimization  technique,  this  kind  of 
sensitivity  analysis  may  be  completed  at  a very  small  cost  in 
computation  time.  It  will  now  be  practical,  therefore,  to  include  the 
results  of  MAXIM-IOS  analyses  of  specific  proposed  arrays  and  array 
patterns  as  a factor  to  be  taken  into  account  during  the  design  process. 

Maxim  compares  well  with  simulation-based  sensitivity  analysis. 
Although  at  first  glance  the  sensitivity  measures  which  can  be  estimated 
through  simulation  seem  to  give  more  information,  these  values  are  in 
fact  dependent  on  a priori  assumptions  concerning  error  distribution. 
MAXIM,  on  the  other  hand,  requires  no  such  assumptions.  The  cost 
factor,  however,  is  the  chief  advantage  of  MAXIM  over  simulation.  This 
advantage  becomes  steadily  more  pronounced  as  array  size  increases. 
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APPENDIX  A 


THE  DIRECTIVITY  FUNCTION  FOR 
A PLANAR  SYMMETRIC  ARRAY 


A. 1 Derivation  From  the  General  Case 

For  a planar  array  In  the  X-Z  plane  consisting  of  n 
magnitude  of  the  contribution  to  the  directivity  pattern, 
angle  of  0,  <),  of  an  element  located  at  is: 

A1exp[j(ui  + c^)]  (A.  1 ) 

where: 

u^  ■ kXj8in©cos0  + kz^cos© 

Ai  is  the  amplitude  shading  coefficient, 
is  the  phase  shift. 

k - 2ir/X 

X is  the  operating  wavelength. 

The  combined  output  of  n elements  is  given  by: 

A1«p[J<u1  + Oj)] 


elements , the 
at  an  aspect 


(A. 2) 
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which  Is  equal  to 

Ai(cos  + j sin  c^) (cos  ^ + j sin  u^ 

Consider  the  contributions  from  four  elements , placed  symmetrically 
about  the  center  of  the  array  at  the  points: 

(.xi,-z1),  and  (x^-z^.  Suppose  further  that  these  elements  are 
excited  identically  (i.e.,have  identical  At  and  coefficients ). This 
result  is: 

Ai(cos  ai  + j sin  ai>  [cos  u.^  + j sin  u2  + cos  u2  + j sin  u ^ 

+ cos^u^)  + j sin(-u^)  + cos(-u2)  + j sin(-u2) ] 

(A. 4) 


where: 

u^  - kxisin©co8()+kz1co8© 
u^  ■ -kx.jSinecoslH-lc^cose 

Since,  for  any  angle  r,  cosr  ■ cos(-r),  and  sin  r “ -sin(-r),  this 

becomes: 


2A  (cos  a1  + j sin  c^)  [cos  ^ + cos  Uj] 

(A.  5) 

By  applying  the  double  angle  formula, equation  A. 5 can  be  stated  as: 

4At  (cos  a±  + j sin  c^)  cos(icxi  sin  6 cos  $)  cos(kz1  cos  0) 

(A. 6) 
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APPENDIX  B 


THE  HOOKE  AND  JEEVES  PATTERN  SEARCH  ALGORITHM 


B»1  Notation 

Before  describing  the  pattern  search  logic.  It  Is  necessary  to 
define  several  parameters  and  notatlonal  conventions. 

1.  ¥ is  a vector  of  n decision  variable  values. 

2.  dj  is  the  perturbation  step  size  for  the  jth  search 
Iteration. 

3.  tjfc  is  the  trial  point  obtained  after  perturbing  about  the 
kCk,  (k«l,n),  decision  variable  in  the  J1"*1  search  iteration. 

4.  & Is  the  search  acceleration  factor. 


B.2  Pattern  Search  Logic 

The  pattern  search  proceeds  from  an  initial  base  point  a,  which 
then  becomes  the  trial  point  for  the  first  Iteration  when  no  variables 
have  yet  been  perturbed,  written  t^.  Each  variable  Is  successively 
perturbed  by  + d^.  t10  and  d^  are  user  specified  parameters.  If  tlk  is 
preferred  over  t^k  t*1®  search  will  move  to  t^k  and  perturb  the 
k*lCh  variable.  The  algorithm  evaluates  trial  points  by  simply 
substituting  the  associated  decision  variable  values  Into  the 
achievement  function.  After  all  n variables  have  been  perturbed,  t^  is 


compared  to  t If  is  not  preferred  over  t^,  the  perturbation 
step  size  is  reduced,  and  the  variables  are  again  perturbed  about  t q. 
If,  however,  t^n  is  preferred  over  t^,  the  search  proceeds  to  a new 
trial  point,  denoted  t g and  given  by: 


C20  “ (tln  " t10) 


i 

t2Q  is  then  compared  to  t^.  If  t2Q  is  preferred  over  tln,  the 
next  iteration  begin*,  from  there.  If  not,  then  t^  becomes  t2g,  d2  is 
reduced  to  a fractional  part  of  dj,  and  the  next  iteration  begins  at 
that  point. 

The  pattern  search  can  terminate  in  the  following  ways: 

1.  A prespecfied  maximum  number  of  search  iterations  may  be 
equalled. 

2.  A user  defined  maximum  number  of  perturbation  step  size 
reductions  may  be  reached. 

3.  The  solution  may  be  within  some  tolerance  factor  of  the 
theoretically  optimal  point. 

4.  The  Improvement  between  two  successive  trial  base  points 
may  be  below  a specified  tolerance  limit. 
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APPENDIX  C 


LINEAR  TECHNIQUES 


r 


C.l  Properties  of  Linear  Models 

A linear  optimization  model  has  the  following  general  form.  An 
objective  function 

m 

I d x 
1-1 

Is  to  be  maximized  subject  to  n constraints: 


i^1ciJ  xi  (>}  bi 


J-l.' 


where  the  xA  are  decision  variables,  and  the  c^ 


are  numerical 

coefficients.  If  a non-linear  optimization  problem  may  be  transformed 
Into  or  approximated  by  a linear  problem.  It  can  be  solved  using  linear 
programming  techniques.  Theoretically,  the  various  linear  programming 
algorithms  always  find  the  optimal  solutlon(s)  If  they  exist.  However, 
in  the  case  of  large  linear  programming  problems  requiring  the 
performance  of  a great  many  operations,  round-off  error  accrued  during 
the  course  of  the  calculations  may  lead  to  an  invalid  solution,  or  cause 
the  algorithm  to  fall  and  find  no  solution  at  all. 

It  is  possible,  by  manipulating  the  directivity  function  and 
aggregating  variables,  to  achieve  a linear  formulation  of  array 
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directional  response.  The  derivation  of  this  result  will  be  reviewed  in 
this  section,  and  the  characteristics  of  the  linear  programming  problems 
so  obtained  will  be  discussed. 


C2.  Linearization  of  the  Directional  Response  Function 

As  shown  in  Appendix  A,  equation  A. 2,  the  combined  output  from  n 
elements  in  a given  direction,  0,9»  is  given  by: 

=i-lAi  6XP[J(Ui  + 

where  uA  depends  on  the  operating  wavelength,  the  element  position,  and 
the  aspect  angle.  For  an  element  in  the  X— Z plane,  for  instance, 

u^  = kx^  sin  0 cos  <|>  + cos  0 


The  output  expression  may  be  equivalently  written: 
n 

Pni  = E (cos  u,  + j sin  u ) A (cos  a + j sin  a ) 

0({l  I XXI 

Expanding  and  rearranging  terms  gives: 
n 

PQ<()  - E [cos  ut  (At  cos  0^)  - sin  ut  (A±  sin  0^)  ] 
n 

+ j E [cos  ut  (A±  sin  c^)  + sin  u±  (A±  cos  0^)] 

If  the  terms  AjCOS  <*  t and  AjSln  ° ± are  understood  to  be  decision 
variables  (i.e.,  the  shading  coefficients  are  now  complex),  and  the 
terms  sin^  and  cosu1  to  be  coefficients,  this  function  assumes  the 
general  linear  form.  The  values  taken  on  by  the  synthetic  decision 
variables,  A^ino^  and  A^os  a ^are  notin  themselves  of  Interest. 
However,  as  long  as  both  of  them  appear  in  the  real  component  of  the 
response  function,  it  will  be  possible  to  recover  the  separate  k±  and  “ ± 
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upon  solution. 

Since  the  coefficients  cos^  and  slnu^  depend  only  on  element 
positions  and  the  6,()  points  at  which  the  pattern  Is  to  be  evaluated, 
numerical  values  may  be  determined  and  substituted  Into  the  response 
function.  Then  the  problem  may  be  solved  using  any  of  a number  of 
computer  packages  for  linear  probrammlng  problems. 


C3.  Linear  Models  for  General  Arrays 

In  1961,  McMahon  et  al,  fl8]  presented  the  transformation  of  the 
directivity  function  Into  a linear  form,  as  reviewed  above,  and 
presented  results  for  several  linear  and  planar  arrays  in  the  12  to  70 
element  range.  This  formulation  applies  to  the  most  general  case  of 
array  design.  The  array  may  be  phased,  and  element  positions  may 
conform  to  any  surface. 

The  chief  disadvantage  of  this  formulation  Is  the  size  of  the 
resultant  linear  programming  problem.  Four  Inequality  constraints  are 
required  at  each  9,<)  point  defining  the  desired  pattern  In  order  to 
create  upper  and  lower  bounds  on  both  the  real  and  imaginary  components 
of  the  response  function.  Furthermore,  4 non-negative  decision 
variables  are  generated  for  each  element  of  the  array.  Thus  the  size  of 
the  resultant  linear  programming  problem  Is  given  by  the  polynomial,  T • 
4n  x 4m,  where  n is  the  number  of  points  needed  to  define  the 
pattern,  and  m is  the  number  of  elements  in  the  array.  Adequate  pattern 
definition  will  usually  require  upwards  of  65  ©,l>  points  per  orthont. 
Thus,  a conservative  estimate  for  the  linear  programming  problem 
generated  by  Juat  a 100-element  (conformal,  planar,  or  linear)  array, 
whose  response  Is  specified  over  a hemisphere.  Is  400  variables  and 


1,040  constraints. 

Unfortunately,  there  is  no  direct  relationship  between  the  sire  of 
a linear  programming  problem  and  the  complexity  function  of  the 
algorithm  required  to  solve  it.  It  has,  however,  been  the  experience  of 
the  author  that  the  use  of  this  method  to  program  arrays  as  small  as  64 
elements  is  prohibitively  expensive. 

The  general  array  linear  formulation  suffers  its  severest 
shortcomings  when  applied  to  arrays  whose  elements  are  symmetrically 
placed.  As  shown  in  Appendix  A for  the  planar  array  case,  such  symmetry 
conditions  may  be  incorporated  into  the  directional  response  function. 
This  type  of  reformulation  allows  the  pattern  to  be  specified  on  a 
fractional  part  of  the  region  over  which  it  is  to  be  defined.  The 
number  of  9,9  points  needed  to  specify  the  pattern  produced  by  an  n 
element  planar  array  displaying  two-way  symmetry  is  half  that  of  a non- 
symmetric  n-element  array.  Four-way  symmetry  halves  this  number  again. 
Since  each  9,9  point  has  4 associated  constraints,  this  method  can 
greatly  reduce  problem  size.  However,  as  shown  in  equation  A. 5,  the 
incorporation  of  symmetry  conditions  causes  the  jsinu^  terms  to  sum  to  0 
and  drop  out  of  the  directivity  function.  The  synthetic  decision 
variables  then  become  indeterminate  in  that  unique  values  for  A^  and  a ^ 
can  no  longer  be  recovered.  Thus,  the  technique  is  Incompatible  with 
the  linear  formulation,  and  symmetric  arrays  will  generate  the  same  size 
problems  as  non-symmetrlc  arrays.  There  is,  however,  one  case  where 
symmetry  can  reduce  problem  size  for  the  linear  formulation.  When  the 
array  is  planar-symmetric  and  relatively  small  (n  < m/2),  the  number  of 
constraints  will  be  reduced  somewhat  if  directional  response  is 
specified  in  only  one  orthont,  and  the  decision  variables  corresponding 
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to  the  elements  in  one  quadrant  of  the  array  are  set  equal  to  those  In 
the  other  three.  For  Instance,  if  the  100-element  array,  260  9,<)  point 
case  given  above  is  treated  in  this  manner,  the  associated  linear 
programming  problem  has  400  variables  and  860  constraints. 

C.4  The  Planar-symmetric,  Non-phased  Case 

Wilson  [20]  has  presented  results  for  planar  arrays  with  two-way 
symmetry,  and  all  element  phase  shift  coefficients  set  to  0.  In  this 
case,  the  response  at  9,1)  from  two  elements  placed  symmetrically  in  the 
X-Z  plane  and  excited  identically  is: 

2A^  cos  u^ 

where  u^  is  defined  as  above.  This  result  can  be  extended  to  the  four- 
way symmetry  case,  and  becomes: 

4A^  cos  u^ 

Then  the  output  of  s sets  of  4 symmetric  elements  is: 

e 

E 4A.  cos  u 
i-1  1 1 

This  function  contains  no  Imaginary  component,  and  since  the  decision 
variables,  A^,  lie  by  definition  between  0 and  1,  only  one  decision 
variable  is  generated  per  set  of  4 elements.  The  100-element  planar- 
symmetric,  non-phased  array  with  the  same  pattern  as  earlier  examples, 
f has  an  associated  linear  programming  problem  of  25  decision  variables, 

and  about  65  constraints.  In  general,  the  size  of  the  problem  in  this 
case  increases  as  the  polynomial,  n/4  x (m/4  + d/4)  where  d is  the 

number  of  9,()  points  in  the  main  beam. 
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The  most  general  objective  function  for  both  McMahon's  and  Wilson's 
linear  programming  problems  minimizes  a weighted  sum  of  deviations  from 
the  specified  responses  at  the  0,  $ points.  Other  objective  funtions 
are  possible,  however.  Wilson's  formulation  [20]  , for  instance, 
minimizes  the  maximum  deviation  over  the  set  of  9,  ^ points. 

The  recent  development  of  a code  to  solve  massive  linear  goal 
programming  problems  [17]  allows  the  reformulation  of  these  linear 
models  as  linear  goal  programming  models.  The  multi-objective  nature  of 
array  synthesis  problems  may  now  be  represented  in  a preemptive  priority 
structure,  as  well  as  in  the  relative  weights  assigned  to  the  goals. 


C.5  MAXIM  for  Planar  Symmetric  Non-Phased  Arrays 

It  is  only  for  the  planar-symmetric  non-phased  array  that  the 
linear  model  has  proved  to  be  of  practical  value,  since  here  the 
associated  linear  programming  problems  are  of  a size  which  is 
economically  feasible.  The  same  holds  for  MAXIM  applications.  The 
linear  version  of  MAXIM  is  as  follows: 

MIN  z'  = E w .d. 

i i i 

n 

St  : E cos  u ■ A + n.  - p.  * b 
^-J  J 1 1 

for  the  i 9,  <)  points  in  the  feature  under  consideration.  d is  a 
negative  deviation  variable  for  the  main  beam,  and  a positive  deviation 
variable  elsewhere.  The  minimization  takes  place  subject  to  a set  of 
constraints  given  by: 

n 


^co.  ut)  + n4  - ?l 


i-l.P 


AJ  ^ Aj  + £j 
Aj  " ^ " £j 


j-1,  n/4 
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where  A*  is  the  design  soltlon, 

and 

2 e 

is 

the  error 

range.  The 

analysis  then  proceeds  as  with 

the 

non- 

linear 

model . 

All  results 

established  for  the  non-linear 

case 

hold 

here 

as  well 

, with  two 

exceptions.  Only  those  error  sources  which  are  resolvable  Into  error  In 
amplitude  settings  may  be  treated  with  this  formulation.  Thus,  for 
instance, error  in  element  placement  could  not  be  analyzed.  Ulth  regard 
to  the  solution  technique  Itself,  however,  the  optimal  solution  Is  now 
guaranteed. 
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